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PREFACE 


The purpose of this book is to provide practical instruction in the 
basic principles of analogue computing by following a course of 
selected problems. The choice of such problems is based on experi¬ 
ence, gained over several years, in teaching students of engineering 
and science in technical colleges and grammar school sixth form 
pupils. 

Successful analogue computation requires a good knowledge of the 
qualitative behaviour of solutions, together with a full appreciation 
of the capabilities of computing units. The theoretical background of 
such knowledge is referred to only briefly as it is felt that it is more 
than adequately treated in the books listed in the bibliography. 
Instead, emphasis has been placed upon the practical techniques 
which need to be mastered in order to make full use of the 
analogue computer. It is hoped such mastery will be attained by 
solving the problems which are arranged in the form of laboratory 
experiments. A routine procedure has been developed to solve these 
problems and, in^ach experiment, sufficient detail is given to enable 
students to programme the computer. At the end of each chapter, a 
supplementary list of problems and answers is given. 

As there is a large number of different types of electronic analogue 
computers in use, no attempt has been made to describe the operation 
of these machines. It is expected that students will become familiar 
with their own machine and have access to the appropriate operator’s 
manual. The reference voltages available in computers depend upon 
whether the machine is a thermionic valve or a transistorized model. 
Accordingly, a standard method has been adopted of using normal¬ 
ized voltages, i.e., expressing all voltages as a fraction of the reference 
voltage. This enables the various methods described in the text to be 
directly applicable to all analogue computers. 

The material in the book covers the needs of most first courses 
in analogue computer programming at college level. Care has 
been taken to restrict the problems to those involving as few com¬ 
puter units as possible, thus obviating the need for large capacity 
computers. It is felt that, if the theoretical side is to be backed by 
practice in the laboratory, then the major part of the course will 
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embody the work on linear operations involved in solving linear 
differential equations with constant coefficients. Usually, in the time 
remaining, only a little work on non-linear operations is possible 
if the student is set the task on the computer himself. The authors, 
therefore, decided to confine their attention solely to the various 
uses of the multiplier and have not discussed the use of further non¬ 
linear devices. 

Since the representation of physical systems by transfer functions 
plays an important role in industrial computation, an elementary 
review of the principles involved has been included in chapter 6 with 
several supporting examples. 

Regarding nomenclature, all resistance values, R, used in diagrams 
are given in megohms, and capacitances, C, in microfarads. Conse¬ 
quently, all products, RC, occurring have the appropriate dimensions, 
i.e., seconds. The symbols used in the figures are shown on page xi 
and, although there is still a variety of so-called ‘standard symbols’, 
the authors have listed those they have found in practice to be the 
most commonly used. 
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Symbols Representing Basic Computer Units 


Component Symbol Bquivalent Operation 



Note: All voltages are in Machine Units, i.e. normalized. The 
electronic multipliers and function generators are available 
in several forms and it is usually sufficient to use the symbols 
given in the second column. 














CHAPTER 1 


Basic Linear Units 

1.1 The Operational Amplifier 

The most important computing unit is the operational amplifier. This 
consists of a d.c. amplifier together with associated input and feed¬ 
back impedances (see Figs. 1.1 and 1.2). The d.c. amplifier is a 
thermionic valve or transistorized unit capable of amplifying voltage 
signals. Several mathematical operations may be performed using 
the operational amplifier. The principal ones are (i) summing a 
number of variable voltages and (ii) producing an integral with 
respect to time of this sum. 

1.2 Use as a Summer 

The d.c. amplifier is usually designed to have a negative gain, — \i 9 
where /i is very large, being approximately 10 8 at zero frequency. If 
e is the voltage at the amplifier input grid, the output voltage e 0 is 

e 0 = -j us. (1.1) 



The current flowing into the grid of the first stage of the amplifier 
shown in Fig. 1.1 is so small that it may be neglected compared with 
the currents i u i l9 z 3 and i f . Hence KirchhofFs first law gives 

+ z*2 + z 3 + if = 0 (1.2) 

and hence by Ohm’s law 

Oi ~ £ ) + (e 2 ~ f) , (fa ~ e) + (go ~ e) _ Q (1 3) 
Ri R 2 Rs R f 
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Inherently e u e 2 , e 3 , e 0 , are voltages of the same order of magni¬ 
tude. As n is very large compared with unity it can be seen from 
Eq. (1.1) that s is very small compared with e 0 , and hence with the 
other input voltages. Assuming that s is zero a very good approxima¬ 
tion to Eq. (1.3) is 


giving 


e l + e 2 ^ e 3 

R i Ri R^ 


Rf 


e 0 



+ b + 


Rt 

R 




(1.4) 

(1.5) 


The quantities R f IRi, R f IR 2 , etc., are referred to as the gains 
associated with the voltages e l9 e 2 , etc., respectively. Some computers 
have a built-in system of resistors providing fixed gains, whilst others 
are provided with external plug-in resistors from which a wide range 
of gains may be obtained. The number of inputs available usually 
varies from about five to ten. 

Attention is drawn to the presence of the minus sign which results 
from the use of an operational amplifier. If the feedback resistor 
and all input resistors are equal a simple sign reversal of the sum 
of the input signals is produced. 


1.3 Use as an Integrator 



In Fig. 1.2 if Q is the charge on the capacitor C at any instant, then 


'V = ^r = cJ r( e ° ~ £ > ^- 6 > 

at at 

and 

h + h + *3 + if = 0, (1.7) 


ej - e e 2 - e 

Ri R2 


e 3 - e 
R 3 


+ cl(e o -e) = 0. 

at 


( 1 . 8 ) 


thus 
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Again assuming that s = 0, 


£L + f* + f*+c^? = 0. (1.9) 

R t R 2 R 3 dt 


Solving for de Q jdt and integrating with respect to t gives the output 
of the integrator at any instant t to be 



(^ + J2_ + 
\CRt CR 2 



dt + e o (0) 


( 1 . 10 ) 


Note again the sign reversal and the gains associated with each 
input voltage, i.e. 1 jCR u 1 ICR 2 and 1/Ci? 3 . The term e o (0) appearing 
in Eq. (1.10) represents the initial value of e 0 . All integrating units 
have provision for inserting the initial value of the output voltage. 
This usually involves feeding the reference voltage via a potentio¬ 
meter, set at the appropriate value, to a specific socket on the inte¬ 
grator. It must be emphasized that the correct initial voltage is 
obtained by measuring the output of the integrator and adjusting 
the potentiometer to give the desired voltage. 


1.4 The Computing Potentiometer 

The computing potentiometer has two uses. First, by applying the 
computer reference voltage to the potentiometer input, any fraction 
of this voltage may be obtained, and second, any variable voltage 
may be multiplied by a positive constant less than or equal to unity. 

Care must always be taken when setting potentiometers to correct 
for loading errors. These become significant when the potentiometer 
load resistance is not large compared with the resistance of the 
potentiometer. 

Thus to perform the operation e t = 0-3^, where e R is the refer¬ 
ence voltage, it is not sufficient to set the potentiometer and obtain 
0-3e R at its output as shown in Fig. 1.3 (a). 



The correct procedure is to set up the complete circuit and adjust 
the potentiometer so that its output under load conditions reads 
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0-3e R . Figure L3 (b) shows an example where the potentiometer is 
followed by an integrator, the point P being where voltage measure¬ 
ment is made. 

In the more modern types of equipment this potentiometer setting 
is easily accomplished using a built-in nulling device. 

Too much care cannot be given to the question of loading, as 
incorrect setting of potentiometers is a major source of error in 
analogue computation. 

The reference voltage in computers is usually ± 100 volts or ± 10 
volts; the former being found in thermionic machines and the latter 
in transistorized equipment. 

To standardize programming the following convention will be 
used throughout the rest of this book. 

The reference voltage will be referred to as 1 Machine Unit 
(M.U.) and all other computing voltages will be taken as fractions 
of this. Thus in a 100-volt machine a computing voltage of 82 volts 
will be 0-82 M.U. and in a 10-volt machine 7*3 volts will be 0*73 M.U. 

1.5 Experiment 1. Use of Summers 

1. Set up an operational amplifier to act as a summer. This is 
usually accomplished by making good a number of connecting links. 
Details are given in the appropriate computer manual. 

2. From the computer reference voltage obtain, with the aid of 
potentiometers, three voltages 0-05, 0T and 0T2 M.U. 

3. By suitable choice of input and feedback resistors generate a 
voltage, 

E = — (lOx + y + z) (1.11) 

where x, y, z are 0*05, 0T, 0T2 M.U. respectively. The circuit is as 
shown in Fig. 1.4. By direct calculation E should be —0*72. 

4. Generate 

E = x — lOy + z. (1.12) 

Equation (1.12) may be written E = — (—x + 10 y — z) from 
which it may be seen that — x, +y, — z are the required input voltages 
(see Eq. (1.5)). The voltages — x and — z may be obtained by supply¬ 
ing the appropriate potentiometers with — 1 M.U. 

5. Generate 

E= — (6x + 1-3 y + 2-4 z). (1.13) 


The circuit is shown in Fig. 1.5. 
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10 10 



When the gain factor required is non-standard, e.g. 1*3, use the 
next higher standard value, 10 say, and precede this input by a 
potentiometer adjusted to offset this increased gain. The setting will 
be 0T3 in this case. The overall gain is the product of the potentio¬ 
meter ‘gain’ and the summer gain. 

Theoretically, potentiometers 1, 2 and 3 will be set at 0-6, 0T3 and 
0*24, but to avoid loading errors it is essential to carry out the test 
procedure described earlier for computing potentiometers. 

6. Generate 

(i) E = — (8x + VIy + 2*lz) (1.14) 

(ii) E = 4-3* - VIy + 5-3z (1.15) 

(iii) E = -lx - 5y + 2-8z. (1.16) 

7. Check the computed values of E by direct calculation. 

1.6 Experiment 2. Use of Integrators 

1. Set up an operational amplifier to act as an integrator. (Consult 
the operator’s manual for the appropriate connections.) 

2. From the computer reference busbar obtain with the aid of 
potentiometers, the three quantities 0-05, 0T and 0-2 M.U. (referred 
to as x, y and z) respectively. 

3. Generate 

E = —J'o-1 dt M.U. (1.17) 

Using the circuit of Fig. 1.6. 



B.A.C.T.— 2 


Fig. 1.6 
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(Note: (i) The initial value of E is zero, i.e. the initial condition of the 
integrator is set to zero, (ii) The input to the integrator, 0*1 M.U., 
should be checked under loading, (iii) In machines where integrators 
have fixed gains, an input provided with a gain of unity is selected.) 

4. By means of the initial conditions facility already mentioned, 
generate 


5. Generate 

(i) E 


E = 10 -J o- 

-T. 

-r. 


1 dt M.U. 


(ii) E 

(iii) E = -J‘ 


(x + y + z) dt. 

(— lOx + y + 0-lz) dt. 
(3x + 0-5y + 0-08z) dt. 


(1.18) 

(1.19) 

(1.20) 
(1.21) 


(iv) 


E = - 0-25 - J* (3-2x - 0-4 y + 0-6z) dt. (1.22) 


6. In each case, apply the signal E to the input of a pen recorder 
calibrated so that 1 M.U. = full-scale deflection. From the known 
values of x, y and z check, by calculation, the computed values of E. 

In the sections to follow full use will be made in diagrams of the 
standard computing symbols shown in the frontispiece. This system 
makes for easier understanding of the more complex computer 
circuits whilst providing all the necessary information. 


1.7 Experiment 3. Simple Applications on Integration 



Fig. 1.7 


1. Using the circuit shown in Fig. 1.7, generate the functions 
y = — 0-U M.U. and y = (M)lt 2 M.U. where 0 < t < 10 s. 

Record the two functions by connecting the voltages representing 
— 0T t and 0-01t 2 to the inputs of a pen recorder calibrated so that 
1 M.U. gives full-scale deflection. 

Some computers have a built-in circuit termed the ‘hold’ circuit 
which can arrest the computation at any predetermined instant. This 
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provides a means of measuring the output signals without incurring 
the dynamic errors of a continuous reading voltmeter. 

The operator’s manual will give full details of the use of the ‘hold’ 
circuit. 

The accuracy of the solution can be checked arithmetically. 

2. A practical application of the above is the examination of the 
motion of a particle travelling with a uniform acceleration. If this 
acceleration is 0*01 ft/s 2 and the initial velocity is 0*05 ft/s the 
velocity v ft/s and the distance s ft after any time t seconds are given 
by 

v = 0*05 + 0*0 It, (1.23) 

s = 0*05/ + £.0*0 U 2 . (1.24) 

As 

v= J* a dt + 0 05 (1.25) 

where a is the acceleration and 



(1.26) 


a circuit arranged as shown in Fig. 1.8 will produce the desired 
solution. 



Note that the output of amplifier 1 is —17 and the solution should 
be interpreted accordingly. 

Representing a velocity of 1 ft/s by 1 M.U. and a distance of 1 ft 
by 1 M.U. record the signals —v and s on a pen-recorder calibrated 
so that 1 M.U. gives full-scale deflection. Once again the hold 
facility may be utilized as an alternative means by which the solution 
may be measured. 

The accuracy of your solutions may be checked by putting values 
of t in Eqs. (1.23) and (1.24). 



CHAPTER 2 


Amplitude Scaling 


2.1 Scale Factors 

The magnitude of each dependent variable of a problem is repre¬ 
sented by a voltage at the output of a computing element. This must 
never exceed 1 M.U. otherwise overloading occurs, warning of which 
is usually indicated either visually or audibly on the computer. 

The constant relating the magnitude of the voltage to one unit 
of the problem variable is called the amplitude scale factor. For 
instance in a problem involving a velocity which varies from 0 to a 
maximum value of 25 ft/s a scale factor K is introduced so that 

1 M.U. = K. 25 ft/s, 
i.e. K = i 5 M.U./ft/s. 

The output of the computing element is labelled v/25 where v is 
the velocity. Any voltage at this output, say 0*6 M.U., is translated 
into the corresponding velocity by multiplying by 25. Hence in the 
given case the velocity is 0*6 x 25 = 15 ft/s. 

The approximate maximum values of the problem variables are 
usually obtained from a knowledge of the physical system under 
examination. The scale factors for each variable are then obtained 
as the ratio 


1 M.U./Estimated Maximum Value. 

In practice these estimates are often too high or too low, a fact 
which is shown by the output voltages from the various computing 
units. It may thus be necessary to re-scale the problem several times 
before a satisfactory solution is obtained. 

These principles are outlined in the following experiments. 

2.2 Experiment 4. A Simple Problem in Dynamics 

A particle moving under a uniform acceleration of 2 ft/s 2 starts with 
a velocity of 10 ft/s. The velocity v and distance s after time t seconds 
are given by 
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v = It + 10 (2.1) 

5 - t 2 + 10/ (2.2) 

where i; is in ft/s and s in ft. 

Obtain graphs of v against t and s against t over a period of 4 s, 
using the computer and suitable recording apparatus. 

1 . Before scaling these equations the maximum values of v and s 
must be estimated. In this simple case they can be calculated exactly 
as 18 ft/s and 56 ft respectively. 

2. The scale factors to produce a maximum of 1 M.U. at the out¬ 
puts representing v and s are given respectively by 

1 M.U. = 7^.18 and 1 M.U. = K 2 .56, 
i.e. K x = iV M.U./ft/s and K 2 ^ ^ M.U./ft. 

These values for K x and K 2 > while permitting full use of the com¬ 
puter dynamic range, would necessitate awkward arithmetic com¬ 
putation in the interpretation of the computer voltages. This may be 
prevented by rounding off the scale factors to slightly lower values. 
Although this does not make full use of the dynamic range the 
resultant loss in accuracy is usually negligible. 

A suitable choice of scale factors is 


K t = 1/20, K 2 = 1/60. 



Fig. 2.1 


3. The computer circuit for this problem is as shown in Fig. 2.1. 
Note that the voltages generated on the machine are labelled v/20 
and —s/6 0. The scaled equations connecting the variables are now 

v/20 - f/10 + 1/2 (2.3) 

s/60 = t 2 /60 + t/6. (2.4) 

These are the Machine Equations and are obtained by multiplying 
Eqs. (2.1) and (2.2) by K x and K 2 respectively. 
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4. As in Experiment 3 the simulation of Eq. (2.3) requires the use 
of an integrator. Its gain G x and initial condition setting I x may be 
calculated as follows. 

As the output has to be u/20 then 



~ G i 


m t 

c x dt + 1 1 
o 


but from Eq. (2.3) 


v 

20 




therefore 


G l e 1 = 


--M.U. 

10 


and 


I, = - M.U. 
2 


(2.5) 


The product G x ei may be chosen in many different ways; in 
Fig. 2.1, e t has been taken as — 1 M.U. and G x as —. 

The distance s is related to the velocity v by the equation 


s 


r. 


v dt + s 0 


where s 0 is the initial value of s (zero in this example). 

As vj20 is available as an input, simulation of Eq. (2.4) requires 
the use of a further integrator with initial condition zero. Its operation 
is described by 



The minus signs take account of the inherent sign reversal of the 
integrator whose gain is seen to be 

Assuming only standard gains of 0T, 1-0, 10-0, etc., to be available 
this value may be obtained by selecting a gain of 1 -0 and preceding 
integrator 2 by a potentiometer set at 0-333 (remember to allow for 
loading of this potentiometer). 

5. The desired graphs will be obtained by recording the outputs 
from integrators 1 and 2 on a chart so calibrated that 1 M.U. gives 
full-scale deflection. This will be equivalent to 20 ft/s and 60 ft for 
v and s respectively. 
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Check these graphs by inserting one or two values for t in Eqs. 
(2.1) and (2.2). 


2.3 Experiment 5. The Discharge of a Capacitor 

A capacitor C has initial charge Q 0 and is discharged from time 
t — 0 through a resistance R. Determine the current flowing in the 
circuit after 6 s in the case where R = 3 M, C — 2 fiF and Q 0 = 100 
coulombs. 

1. The laws relating the voltage V, current /, and charge Q at time 


t seconds are 

i = = VjR 

dt 

(2.7) 

and 



V = QIC. 

(2.8) 

Hence 

- = - i/RC. 
dt 

(2.9) 

Initially, 



h — QolRC- 

(2.10) 


2. Since energy is lost the maximum numerical value of the current 
is 

Qo/RC = 100/6 amp. 


The maximum scale factor for i is therefore 


K = M.U./amp. 


To simplify circuit details and to facilitate graph scales take 
K = • 

3. The scaled equation corresponding to Eq. (2.9) is thus 



Solution of Eq. (2.11) on the computer may be obtained using the 
circuit of Fig. 2.2. 


This is obtained by first considering the output of the integrator. 
If this is taken as if 20 then feeding back this voltage and integrating 
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according to Eq. (2.11) the desired result will be achieved. To 
accomplish this, a gain of 1-0 is used in the integrator in conjunction 
with a potentiometer set at 0T67. The initial value of i being 100/6, 
that of z/20 is 5/6 and is obtained using a potentiometer set at 0*833. 



4. Record the solution over a period of about 10 s on a recorder 
chart calibrated to give a full-scale deflection for 1 M.U. The value 
of the current flowing after 6 s can then be read from the graph. 

5. Obtain the solution to the same problem with C half its original 
value. (Hint: Choose a new scale factor for i and adjust the potentio¬ 
meters accordingly.) 

2.4 Experiment 6. A Problem in Particle Dynamics with 
Variable Acceleration 

A particle moves in a straight line with an acceleration proportional 
to time t and directed away from a fixed point O on the line. The 
acceleration is 1 ft/s 2 when t = Is and the particle starts from a 
point P, such that OP = 25 ft, with a velocity 18 ft/s towards O. 
Obtain the graphs of velocity v and distance x against time t. 
Find v and x when t is 8 s, and the distance covered when the velo¬ 
city is zero. After how many seconds does the particle reach 01 
1. The equations representing the motion are 


II 

"§ l*« 

(2.13) 

" d - V dt - 18 

O 

(2.14) 

j v dt + 25 

(2.15) 

Jo 



where x, v and dvjdt are measured positively in the direction OP. 
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2. In this problem the greatest numerical values of x, v and dvjdt 
can be obtained by direct solution of the equations. If t does not 
exceed 10 s these maximum values are 47 ft, 32 ft/s and 10 ft/s 2 
respectively. 

3. The maximum scale factors are: 


« dv 
for —, 
dt 

= Vo M.U./ft/s 2 

for v, 

K 2 = V-2 M.U./ft/s 

for x, 

K 3 = Vt M.U./ft. 


Suitable working values are: 

K t = Vo M.U./ft/s 2 
K 2 = To M.U./ft/s 
K 3 = Vo- M.U./ft. 

4. The machine equations are: 


1 (Jd 1 

To ~dt ~ 10 


(2.16) 

1 _ _ 1 


18 

(2.17) 

50 17 f 

>Jo 10 dt 

50 

1 

['-±vdt + 

25 

(2.18) 

50 X _ 

Jo 50 

50 


5. The computer circuit to solve these equations is as shown in 
Fig. 2.3. Notice the use of a sign reverser to obtain -f x/50. 

6. Record the outputs from integrators 2 and 4 and hence deter¬ 
mine the solutions of the problem. 


+ 1 MJ 


U. 

^~ 10 W 




+ 1 M.U. 



CHAPTER 3 


The Solution of Problems Involving Ordinary 
Differential Equations with Constant Coefficients 


3.1 Introduction 

Very many physical problems may be represented by mathematical 
models in the form of differential equations. It must be pointed out 
however, that in practice, only approximate representation of 
physical systems by such models is possible. Nevertheless, much 
important information may be obtained from the use of differential 
equations and their solution. 

In order to facilitate greater understanding of the methods of 
solution involved in the experiments of this chapter a brief review 
of the nature of the solutions and method of programming for first 
and second order differential equations follows. 

3.2 First Order Differential Equations 

The general first order differential equation with constant co¬ 
efficients is 

ax + bx = f{t) (3.1) 

where a and b are constants and f{t) is an arbitrary function. 

This equation has the complete solution 

x = Complementary Function (C.F.) + Particular Integral (P.I.). 

The particular integral in general takes the form of f{t). The com¬ 
plementary function is the general solution of Eq. (3.1) when/(t) = 0, 
and is 

x = Ae~l\ (3.2) 

If bja is positive the magnitude of the C.F. decreases with in¬ 
creasing t , while if bja is negative it increases. 

Equation (3.1) may be re-arranged 

X = -~x + - f(t ) 
a a 


(3.3) 
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Hence 

-x - l/(ol dt + x o (3.4) 

a a J 

where x 0 is the initial value of x. 

The computer circuit (Fig. 3.1) used to solve Eq. (3.4) will thus 
consist of an integrator supplied with the appropriate initial condition 
and fed with inputs x and — f(t ). The latter is obtained from a suit¬ 
able circuit (not shown) and the former from the output of the in¬ 
tegrator itself. 




Fig. 3.1 


3.3 Second Order Differential Equations 

The general linear second order differential equation with constant 
coefficients is of the form 

ax + bx + cx = F(t). (3.5) 

When the coefficients a, b, c are all positive this equation represents, 
ideally, many common physical phenomena, e.g. a spring oscillating 
in a damping fluid. 

Equation (3.5) may be written more conveniently in the form 

x 4- 2 L,a) n x + co 2 n x — fit) (3.6) 

where £ and co n are constants which play a significant part in the 
nature of the solution of the equation. Once again the complete 
solution consists of the sum of a particular integral and a com¬ 
plementary function. 

The particular integral, in general, assumes the form of fit), and 
the complementary function is the general solution of Eq. (3.6) 
when f{t) = 0 and is 

(i) x = A sin cot for £ = 0. 

(ii) x = e~^ nt \_A cos oo ny J{ 1 - £ 2 )f + B sin a> n% /(l - £ 2 )t] 

for £ < 1. 

(iii) x = e~ a>nt (At + B) for £ = 1. 

(iv) x = + Be ~ aWtt2 “ 1)f ] for £ > 1. 



16 BASIC ANALOGUE COMPUTER TECHNIQUES 

These results are illustrated in Fig. 3.2, the wave forms being 
termed (i) undamped, (ii) under-damped, (iii) critically-damped and 
(iv) over-damped. 

C which clearly controls the rate of decay of stable solutions is 
called the Damping Ratio. Negative values of C lead to exponential 
terms whose magnitude increase indefinitely with t and hence con¬ 
tribute to unstable solutions. 



co n is the angular frequency in radians per second at which oscilla¬ 
tions occur when C = 0 and it is defined as the undamped natural 
frequency. In the case where 0 < £ < 1 the actual frequency of 
oscillation co is given by 

cu = <0*7(1 - C 2 ). (3.7) 

The actual period of oscillation is 2 n/co seconds. 

Equations such as (3.5) are solved on the computer by successive 
integration as follows. Solving for the highest derivative present 


b. 

x = —~x 
a 

After one integration 


c 

-x 

a 



x = 



where x 0 is the initial value of x. 
Integrating again 




dt + x 0 


(3.8) 


(3.9) 


x = 


x dt + x 0 . 


(3.10) 
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The computer circuit to solve Eq. (3.5) is developed by first 
considering Eq, (3.9). As illustrated in Fig. 3.3 this is accomplished 



Fig. 3.3 


by summing the three quantities appearing on the right-hand side of 
(3.9) and integrating this sum. 

The sign change in passing through the integrator gives +x as 
the output, x is produced from this by use of another integrator. 

The output of the first integrator, being + x, is directly available, 
after suitable multiplication, for the input signal. The output from 
the second integrator must be sign inverted before it can be made 
available for the input -\-(cja)x and the other signal —F(t)/a is 
obtained from within the computer. 

Figure 3.4 shows the completed circuit in unsealed form. 



Fig. 3.4 


3.4 Estimation of Scale Factors 

It must be emphasized that there is no simple method for estimat¬ 
ing the maximum values of the variables and their derivatives in 
differential equations. Consideration must be given both to the 
initial conditions and the value of the function f(t) occurring on the 
right-hand side of the equation. 

For instance in the first order Eq. (3.1) if f(t ) — 0 
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* = --X (3.11) 

a 

and the complete solution is 

_b 

x = x 0 e * (3.12) 

where x 0 is the initial value of x. 

If b\a is positive, the greatest numerical value of x, |x| m , is |x 0 | and 


XL = 


(3.13) 


A more detailed application is given in Experiment 7. 

In the second order Eq. (3.6), when f(t) = 0, C = 0 and x = x 0 , 
x = 0 at t = 0 then 

|*|m = (3.14) 

|*|- = (o 2 „\x\ m (3.15) 

where \x\ m = \x 0 \. 

When 0 < £ < 1 the values of |x| m and |x| m will be less than those 
given in Eqs. (3.14) and (3.15), the corresponding relations being 

|*|m < H*|m (3-16) 

\x\ m < co 2 \x\ m (3.17) 

where co = oo n ^J( 1 — C 2 )- 

As will be shown in subsequent experiments these results provide a 
rule of thumb method for estimating the maximum values of the 
problem variables. 

For a differential equation of any order the following empirical 
rule suggested by A. S. Jackson in his book Analog Computation 
may be used to obtain estimates of the maximum values of x and its 
derivatives. 


3.5 The Equal Coefficient Rule 

Consider the nth order differential equation 

a n x M + a„- 1 x (n_1) + . .. + a y x w + a 0 x =/(<) (3.18) 

where x (r) s — for r = 1,2. n. 

df 

For simplicity take zero initial conditions and 
f{t) = 0 when t < 0 
fit) = A, a constant, when t ^ 0. 
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With these conditions the following results may be obtained from 
Eq. (3.18) 

x** = A/a n (3.19) 


and 


x m = 2Aja 0 


* 


(3.20) 


where x m and xj^ 0 are maxima again referring to absolute magnitude. 

Dividing each derivative by its maximum value, Eq. (3.18) 
becomes on re-writing, 



+ x 


(n-l) 

m 


a n- 




+ 



+ a A x 


U) 

m * 



+ 2 A 


x 

2 A/a 


= /( 0 - 


(3.21) 


This is called the Normalized Equation and the quantities in 
brackets the Normalized Variables. For instance x (1) /x^ I 1) is the 
normalized first derivative and will not exceed the value 1. The 
estimates for the maximum values x^\ x ( ^\ ..., xj£ _1) are then 
determined by making the coefficients of the normalized variables 
equal, with one exception; the coefficient of (x/x m ) will be twice that 
of the other coefficients. 

One restrictive condition applies. The values of x { J^\ ..., x^ n) 
must form an increasing or a decreasing sequence. 

The application of this rule is shown in the following example. 

Consider the equation 

x + 2x + 9x + 18* = fit) (3.22) 

with x = x = x = 0 at t — 0 

and where f(t) = 0 for t < 0,f(t) = 36 for t ^ 0. 

(a) From Eqs. (3.8) and (3.9), x m = 36 


x m = 72/18. 

( b ) The normalized equation thus becomes 



(c) For equal coefficients of the normalized variables (with the 
exception of x/x m ) 

36 = 2x m — 9x m 

hence x m = 4, x m = 4, x m = 18 and x m — 36. 
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These agree very well with the actual maximum values which are: 
x = 3*02, x = 3*85, x = 10*34 and x = 36. 

Once the maximum values of the dependent variable and its 
derivatives have been found, the corresponding scale factors can be 
evaluated and the original differential equation written in a scaled 
form. 

To do this, each problem variable is multiplied by its scale factor, 
and the original differential equation re-written in terms of these new 
variables. Since the scale factor for any variable is calculated from 
the relation 


Scale factor = 1/Estimated Maximum Value 

r \ 

it will be noted that t^ese ne\^ variables are the normalized variables. 
Hence the seated equation ^-equivalent to the normalized equation. 

The scal^equati^j^solved for the highest scaled derivative as in 
Eq. (3.8), oQe_iq|^ration performed, and the L.H.S. rescaled for the 
lower‘ordesJSS^ative. This new equation, suitable for execution on 
the computer^ is the Machine Equation. In the case of the previous 
example (Eq. (3.22)) the scaled equation is 



and the machine equation is 



where (x\ 18) 0 is the value of ic/18 when t = 0. 

In general since the highest order derivative (in this case x) is not 
present in the machine equation, its scale factor is not needed for the 
computer circuit. Only occasionally is a graph of this function against 
t required; in this instance the scale factor is used to enable the 
appropriate function to be generated. 

The following experiments are suggested in order to develop these 
methods of preparing differential equations for solution on the 
computer. 

3.6 Experiment 7. Solution of a First Order Differential Equation 

Given that 


2x + x = 0 


(3.26) 
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and x = 50 when t = 0 obtain a graph of x against t over a period 
of 10 s. 

1. As the coefficient of x is positive the solution of Eq. (3.26) is an 
exponential decay (compare Eq. (3.2)). The greatest numerical value 
of x thus occurs when t = 0, i.e. |x| m = 50. 

Using the relationship given in Eq. (3.13), |x| m — 25. 

2. The scale factors for x and x are respectively 


K 1 — sq and K 2 — 2 1 5 * 


The scaled equation is 


i.e. 


2 x 25| — | + 50(—I = 0, 
\25/ \50J 


— + 
25 


X 

50 


= 0 . 


(3.27) 


3. Solving for x/25 and integrating this gives 


x 

25 


o 50 


dt + 



(3.28) 


But, since it is desired to generate a signal in the computer equivalent 
to x/50 the machine equation is 


x 

50 



(3.29) 


4. The simulation of Eq. (3.29) requires one integrator whose 
output is x/50 and input x/50. If its gain is G t and initial condition I u 
then for this integrator, output 



Comparing Eqs. (3.29) and (3.30) leads to the values 

G, = 1 and /. = (—) =1 M.U. 

2 V50) o 


G x may be obtained by using an integrator gain of unity preceded by 
a potentiometer set to 0-5. 

The output from this integrator may be fed directly into the input 
to the potentiometer, since it is the required function, x/50. Figure 
3.5 shows a computer circuit for Eq. (3.29). 

B.A.C.T.— 3 
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5. Record the output from the integrator and check the accuracy 
of solution using the values in Table 3.1. 

TABLE 3.1 


t 

0 

1 

2 

10 

X 

50 

30-3 

18-4 

0-34 


6. Using the same initial condition x 0 = 50 write out the machine 
equation for 

2x + x = -20 (3.31) 

and modify the circuit of Fig. 3.5 to compute the new solution. 

The analytical solution in this case is 

x = 10e~ t/2 - 20, (3.32) 

from which the graph of x against t may be checked. 


3.7 Experiment 8. Solution of a Second Order Differential Equation 
If 


x + 2x + llOx = 0 


(3.33) 


and x = 25, x = 0 when t = 0, obtain the graphs of x against t and 
x against t ; given that — 25 < x < 25. 

1. The natural undamped frequency of the system represented 
by Eq. (3.33) is 

co n = ^/llO = 1049 rad/s. 


The damping ratio is 

£ = _ 2 — 

2^110 


0-0954. 


The actual frequency of oscillation is 


co = 7110^/(1 - 1/110) = 10-44 rad/s. 
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2. Using the relations (3.16) and (3.17) which are applicable here, 

|x| m = 10-44 x 25 = 261 
and |x| m = 2725 

3. The scale factors are: 


for x, Kl = Ts 
for x, K 2 = t 
for x, K 2 = 2 i 2 s' 


It will be more convenient to take K 2 = 1/500 and K 2 = 1/5000 as 
these figures make for easier interpretation of chart recorded values. 

4. The scaled equation is 

500of JL) + looo /+ 2750j— J = 0. (3.34) 

\5ooo; V 500 / V 25 / 

Solving for x and integrating once 

500(/—) = - f T1000/—) + 2750| — )"1 dt. (3.35) 

Vsoooy JoL \50o; V 25 /J 

(The initial value of x is zero.) 

5. Since it is required to generate a signal equivalent to x/500 
in the computer the machine equation is 


x 

500 



(3.36) 


From this equation it can be seen that there is no need to scale the x 
term as no corresponding signal appears in the computer. 

6. The computer circuit is constructed as described earlier in the 
chapter for Eq. (3.5). Any confusion regarding sign inversion is 
avoided by writing Eq. (3.36) with the negative sign outside the 
integral. Thus to obtain an output of x/500 from the first integrator 
inputs of x/500 and x/25 are necessary. If the respective gain factors 
are G x and G 2 and the initial condition is I u then, 

dt + I x . 

Comparing Eqs. (3.36) and (3.37) gives the values 




= 2, G 2 = 5-5, and /, = 0. 


(3.37) 
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The signal -x/25 is obtained using a further integrator with gain (? 3 
and initial condition / 2 . Then 


Now since 


=- r oidt 

25 Jo \5Q0J 


+ f 2* 


S _ ±\dt~ 1 

25 Jo 125 


(3.38) 

(3.39) 


the values for (7 3 and / 2 are found to be 


G 3 = 20 and I 2 = — 1. 

A sign inverter of gain unity provides the necessary signal + x/25 
for the input to the first integrator. The signal x/500 is available to be 
used as the other input. Figure 3.6 shows the resulting circuit when 
standard gains in association with computing attenuators are used. 



7. Obtain the graphs of x and x against t from the signals at the 
output of amplifiers 1 and 3. 

8. Check the accuracy of your solution from the following in¬ 
formation: 

x: Frequency 10-44 rad/s 

Initial value 25 units 

Final value 0 units. 


TABLE 3.2 


Peak 

1st 

2nd 

3rd 

4th 

5th 

x value 

+ 25 

-18-5 

+ 13-7 

-101 

+ 7-5 

t value 

0 

0-3 

0-6 

0-9 

1-2 
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x: Frequency 10*44 rad/s 
Initial value 0 units/s 

Final value 0 units/s. 


' TABLE 3.3 


Peak 

1st 

2nd 

3rd 

4th 

x value 

-221 

+ 167 

-124 

+ 92 

t value 

0*15 

0*45 

0*75 

1-05 


3.8 Experiment 9. Generation of Sine and Cosine Functions 

Generate 5 sin 21 and 10 cos 21. 

1. x = a sincoMs the solution of the differential equation 

x + co 2 x = 0 (3.40) 

where x = 0 and x = aco at t = 0. 

When co = 2, a = 5 Eq. (3.40) becomes 

x + 4x = 0 (3.41) 

where x = 0 and x = 10 at t = 0. 

2. From Eq. (3.41) the actual frequency of oscillation, which is 
equal to the natural undamped frequency, is 2 rad/s. 

3. Since x = 5 sin2f 

|*|« = 5, |x| M = 10 and |x| m = 20. 

4. The scale factors are 

for x, K 1 = | 
for x, K 2 = T l o 
for x, ^ 

5. The scaled equation is 

20 (i) + "(s) = °' <3 ' 42) 

6. Solving for x and integrating once 


(3.43) 
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7. As it is required to generate a signal i/10 in the computer the 
machine equation is 



(3.44) 


8. Programming the computer from this equation, integrator 1 
has an input of x/5 and an output of i/10. If the gain factor is G t 
and the initial condition I u then 


+ <3 ' 45) 

Comparing Eqs. (3.44) and (3.45) 



The signal —x/5 is obtained using a further integrator with gain G 2 
and initial condition I 2 . Then, since 

Hh-R* (146) 

D J o 3 

“ d -5““£ C! (ro)‘" + /2 (3 ' 47) 

it follows that E 2 = 2, I 2 = 0. 

A sign inverter of gain unity provides a signal x/5 for the first in¬ 
tegrator and the resulting circuit is as shown in Fig. 3.7. 



9. The results obtained may be easily checked by measuring the 
amplitude and frequency of the solution x. 

10. From the same circuit derive and record a signal 10 cos 2 (. 
The functions a sin cot and a cos cot are frequently used in comput- 
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ing, particularly as driving functions. It is recommended that students 
make themselves absolutely familiar with the above procedure. 

3.9 Experiment 10. A Second Order Differential Equation 
with Variable Damping Coefficient 

The displacement x of a body acted upon by certain forces satisfies 
the equation 

x + rx + 9x = 0 (3.48) 

and x = 100 cm, x = 0 when t = 0. r is a variable damping coeffi¬ 
cient, which is a measure of the resistive force suffered by the body. 
Obtain the graphs of x against t for the cases, (i) r = 1-0, (ii) r = 6*0, 
(iii) r = 10-0. 

1. All three cases may be dealt with using a single circuit provided 
a potentiometer is included in the damping loop and made propor¬ 
tional to r. The amplitude scale factors for the variables can be 
calculated using the lowest value of r, since this will, in general, lead 
to the most oscillatory condition and hence the largest amplitudes 
of x, x and x. 

Since r is positive the system is a decaying one and it will be a good 
enough approximation to take 100 as the maximum value of x, 
i.e. x m = 100. 

The natural undamped frequency of the equation is co„ = 3. 

The damping ratio f = r/6 = 1/6 in the case of least damping. 
In this case the actual frequency of oscillation is 

co - 3v(l - ■&) = 2*96 rad/s. 

= 3 approximately. 

2. Using the relations (3.16) and (3.17) 

|x| w = 3 x 100 = 300 cm/s 
|x| m = 9 x 100 = 900 cm/s 2 

3. The scale factors are 

for x, K t = 
for x, K 2 = no 
for x, K 3 = g^Q. 

4. The scaled equation is 



(3.49) 
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5. Integrating this gives 



6. The computer circuit is constructed as described for Experi¬ 
ment 8. Integrator 1 has inputs of x/300 and x/100 and an output 
x/300. if the gain factors on the respective inputs are G u G 2 and the 
initial condition is I u then 


x 

300 


~JI c ' 


300 


+ G : 


100 


dt + / 1# 


Comparing coefficients in Eqs. (3.50) and (3.51) 


(3.51) 


G ! = r, G 2 = 3, I t = 0. 


To accommodate the largest value of r, i.e. 10, a potentiometer 
whose setting is r/10 is used in conjunction with an integrator gain 
of 10 to produce gain G l .G 2 is similarly obtained by a potentiometer 
set at 0*3 followed by an integrator gain of 10. 

The signal — x/100 is obtained by feeding x/300 through integrator 
2. If G 3 and / 3 are the gain and initial condition of this integrator 



A sign inverter provides the signal + x/100 for the input to in¬ 
tegrator L The output of integrator 1 being x/300 is immediately 
available as the other input. Figure 3.8 shows the resulting scaled 
computer circuit. 



Fig. 3.8 
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This problem involves the use of three potentiometers and two 
integrators but in more complex problems it is more usual to employ 
many more of these units. It is therefore essential that a methodical 
check procedure be carried out in the computer before attempting 
to obtain a result. One such procedure will now be outlined and 
applied to this experiment. 

First a static check consisting of two parts is applied. Part one, 
a pencil and paper check, is made before carrying out any patching. 
The procedure is as follows: 

{a) Represent the problem variables and their derivatives by scaled 
initial condition voltages at the output of each integrator. These 
voltages are quite arbitrary but are chosen to produce reasonable 
outputs, they need have no physical significance. 

(b) On the computer diagram write down these selected scaled 
voltages next to the outputs of the appropriate integrators. 

(c) Using these voltages perform the various operations of 
summing, multiplication by a constant, etc., and hence calculate the 
voltages at all the input and output terminations of all the computing 
components. As integrators are not operational in the initial condi¬ 
tion position this calculation is not carried beyond the input stage 
of any integrator. 

id) Substitute the assumed initial condition values of the variables 
into the original (i.e. the unsealed) equation describing the problem 
and calculate the value of the highest derivative. This value should 
tally with the voltage level obtained in (c) for that point. If the 
actual voltage levels are too low or zero it is necessary to change 
some of the assumed values of the variables and repeat the above 
calculations. 

The problem is now patched on the computer and the second part 
of the static check is carried out. This consists of setting the potentio¬ 
meter and the initial conditions to the static check values and then 
measuring the outputs of all the computing components with the 
switch in the Initial Condition position. 

To check that the input resistors to an integrator are of the correct 
value the following procedure may be followed. 

In the given circuit (Fig. 3.9) the gains through the integrator are 
1, 10 and 1. Disconnect junction J from the integrator and connect 
to an amplifier with a feedback resistor of 1. The output of this 
amplifier should be 1(1) + 10(0-02) + (0-1) = 1*3. 

As this produces overloading the gains need to be reduced by a 
factor of 10, i.e. a feedback resistor of 0-1. 
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Having completed the static checks it is often possible to perform 
various dynamic checks; these however depend upon the information 
available and it is obviously impossible to generalize. 



This check procedure will now be applied to Experiment 10. For 
clarity some form of tabulation is recommended (see Table 3.4). 

(1) The potentiometer settings are calculated from the computer 
diagram Fig. 3.10 and entered on the sheet. 

(2) Assume the following values for the variables 

x = 300, x = 100 and let r = 1. 

(3) On a copy of the computer diagram mark in the static check 
voltages. (These appear in the small rectangles on the diagram.) 
Notice that integrator 1 is now supplied with an initial condition for 
the purpose of the check procedure. 



(4) Substitution in the original equation of the problem gives 

x = ~(300r + 900) ft/s 2 . (3.54) 

The input to integrator 1 should thus be 

- — = +(r + 3) M.U. = +4 (3.55) 

300 


in the case r — 1. 
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From the circuit above the input at the grid of integrator 1 is 
10(0T)(+1) + 10(0-3)(+l) = 4. This checks with the calculated 
value and indicates that the original equation is correctly modelled 
by the computer circuit. 

(5) As previously stated the computer is now patched and with 
the above settings the output voltages are checked. 


TABLE 3.4 
Problem Sheet 


Potentiometer 1 

2 

3 

4 

5 

Setting (r) 

0-1 

0-3 

0-3 

1 

1 

Amplifier 

Output 


Static check 




Calculated 

Measured 



Check point 
(M.U.) 

Output 

(M.U.) 

Check point 
(M.U.) 

Output 

(M.U.) 

1 

X 

300 

r + 3 

+ 1 



2 

X 

"100 

3 

~1 



3 

X 

Too 


+ 1 




3.10 Experiment 11. Solution of Simultaneous Differential Equations 
when Knowledge of the Problem Variables is Available 

Solve using the analogue computer the equations 

x — 2y = 20 
y + 2x = 0 

given that x = 50 and y = 40 when t = 0. 

Assume that — 70 < x < 70 and — 80 < y < 60, the analytical 
solution of Eqs. (3.56) and (3.57) being 

y = 50^/2 sin (2 1 - tc/ 4) - 10 
x = - 50^2 cos (It — tt/4). 


(3.56) 

(3.57) 
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(1) Using the analytical solution the maximum values of x and y 
can be calculated as 


|*|. = 100V2, \y\ m = I00yj2. 

(2) The scale factors are 

for x, K x = 1/70. Use 1/100 

for jc, K 2 = 1/100V2. Use 1/150 

for y, K 3 = 1/80. Use 1/100 

for y, A 4 = 1/100/2. Use 1/150. 


(3) The scaled equations are 



and 


Thus 



+ 2001 


x 

Vioo 


= 20 

= 0. 


and 


x 

150 


4/^ 2 

3\100 J 15 


y _ 4/ x \ 

150 ~ ~3^l00y 


(3.58) 

(3.59) 

(3.60) 

(3.61) 


(4) Integrating each equation once and rearranging gives the 
machine equations 



Note again that a minus sign is placed outside each integral. As 
previously stated this is to be compared with the inherent sign 
reversal occurring in the electronic integrators. The signs attached 
to each term of the integrands of Eqs. (3.62) and (3.63) become the 
signs required for the corresponding input connection of the com¬ 
puter circuit. 

(5) Equations (3.62) and (3.63) can be solved using two integrators. 
The first supplied with inputs — y/100 and —1/5 M.U. develops an 
output of x /100 while the second, supplied with an input of x/ 100 , 
produces an output of y/ 100 . 
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If the gain factors on the inputs of the first integrator are G t and 
G 2 and its initial condition is I l9 then its output is 


x 

ioo 


-P ' -\g x (-?L\ + gJ^\ 

Jo L V 100 / \5/_ 


dt + /j. 


(3.64) 


Similarly if <j 3 and I 2 are the gain factor on the second integrator 
and its initial condition respectively then 


y 

ioo 



(3.65) 


Comparing Eqs. (3.62) and (3.63) with (3.64) and (3.65) gives 


h 


Gi 



“ 2, G 2 = 1, 
0-5 M.U., I 2 = 



0-4 M.U. 


A sign inverter of gain unity converts the output of integrator 2 to 
the appropriate sign for the input of integrator I. The complete 
circuit is as shown in Fig. 3.11. 



Potentiometer 1 2 3 4 5 

Setting 0*2 0*2 0*5 0*2 0*4 

(6) Perform the static check procedure described earlier to confirm 
the voltage levels occurring in the initial set condition. 

(7) Record the values of x/100 and j^/100, and test the accuracy of 
solution using the analytical form given above. The frequency of x 
and y is 2 rad/s, thus if a number of complete periods are measured 
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from the graphs of x and y against t the period and hence the angular 
frequency of oscillation may be checked. Do not simply use large 
values of t to test the solution, since small percentage errors in the 
timing can cause large errors in the value of x and y, caused by the 
resulting phase shift of the curve. 

3.11 Experiment 12. Solution of Simultaneous Differential Equations 
when Prior Knowledge of the Variables is not Available 

Obtain the graphical solution of the simultaneous differential 
equations 


dx 

— + 2x + y = 0 
dt 

(3.66) 

Q + x + 2y = 0 
dt 

(3.67) 


given that x = 1 ,^ = 0 when t = 0 . 

(1) In this problem nothing is known of the maximum values 
attained by the variables. However estimates may be made in the 
following way. 

Consider the first equation with y regarded as constant at its initial 
value, i.e. y = 0. The equation reduces to 

— + 2x m 0 (3.68) 

dt 

and has the solution x = e~ 2t using the value x = 1 at t = 0 . 
Consequently the greatest value of x is 1. 

If this maximum value of x is used in Eq. (3.66) in place of the 
variable x, then 

— + 2y = — 1 (3.69) 

dt 

which has solution y = i(e~ 2t — 1 ). 

Hence \y\ m = 

(2) The scale factors are 

for x, K x = 1 
for x, K 2 = \ 
for y, K 3 = 2 
for y, K a = 1 

using the principle outlined in Experiment 7. 
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(3) The scaled equations are: 
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2^) + 2(x) + 1(2 y) = 0 

(3.70) 

\ / 

and 

y + x + 2y = 0. 

(4) The machine equations are 

(3.71) 

X= h 

2x + -(2m dt + 1 

^ J 

(3.72) 

and 



2 y = - 

[2x + 2(2y)] dt. 

0 

(3.73) 


Referring to previous experiments it can be seen that the gain factors 
associated with each amplifier input are the actual coefficients of the 
scaled variables appearing in the machine equation. Thus in Eq. 
(3.72) the x input to the integrator has an associated gain of 2, this 
being achieved using a gain of 10 together with a series potentio¬ 
meter set to 0-2 (potentiometer 1 of Fig. 3.12). In this and future 



Fig. 3.12 


experiments this method will be adopted in programming the 
computer. 

(5) Equations (3.72) and (3.73) are solved in the computer using 
a circuit containing two integrators as in Fig. 3.12. To obtain the 
desired gains the potentiometer settings are as follows 

Potentiometer 12 3 4 

Setting 0-2 0*5 0*2 0*2 
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(6) Perform the static check using some non-zero value of (2y) 0 . 

(7) Record the solutions of x against t and 2 y against t. 

(8) Check the accuracy of solution by direct calculation from the 
solution 

x = \(e~ % + e ~ 3t ) 

y = l(e~ 3 ' - e~‘). 

3.12 Experiment 13. Forced Mechanical Oscillations 

A body of mass 8 lb is hung by a spring producing in the spring an 
extension of 3-2 ft. The upper end of the spring is made to execute 
a vertical oscillation z = 12 sin At , z being measured vertically 
downwards in feet. The mass is subjected to a frictional resistance 
whose magnitude in pounds weight is one-half of its velocity in ft/s. 
Obtain a graph of the displacement x from the equilibrium position. 
Take g to be 32 ft/s 2 . 

It is first necessary to obtain the equation describing the physical 
system. 



Fig. 3.13 


For a body of mass m lb oscillating as stated the forces upon it are 

(i) The gravitational force mg. 

(ii) The spring force which opposes the displacement and is 
proportional to the extension of the spring. 

(iii) The damping force which opposes the velocity and is pro¬ 
portional to the velocity. 

(iv) The excitation z. 

Newton’s second law of motion states that the rate of change of 
momentum is proportional to the impressed force and takes place 
in the direction of the straight line in which the force acts. For a 
constant mass the rate of change of momentum is equal to the pro¬ 
duct of the mass and acceleration. In the given system the equation 
of motion is 






i.e 
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my = algebraic sum of the forces in the y direction, 

my = -Kg(y - z - 1) - rgy + mg. (3.74) 

In equilibrium (Fig. 3.13 (a)) 

0 = —Kge -f mg 

i.e. K = —. 


Substituting for K and rearranging, Eq. (3.74) becomes 
m m 

—y + —(y - 1 - e - z) + r$ = 0. (3.75) 

9 e 

From Fig. 3.13(Z>) it can be seen that the displacement x from the 
initial equilibrium position is y — 1 — e, i.e. 


y - x + 1 + e. 

Substituting for y and its derivatives Eq. (3.75) becomes 


m 

—x + rx + 
9 


m 

-x = 
e 


m 

—z. 

e 


(3.76) 


Substituting the problem values the describing equation is therefore 
x + 2x + lOx = 120 sin At. (3.77) 

1. The maximum values of the problem variables have now to be 
found. In practice a working knowledge of the type of problem will 
in most cases enable a good estimate to be made. For this particular 
problem a value for the maximum values of x and x may be found 
using the equal coefficient rule and treating the right-hand side of 
Eq. (3.77) as a constant equal to 120. 


2 170 

This rule gives |x| m = = 24 (see (3.20)) 

and |x| m = ip = 120. 


The normalized equation is 
+ 2.x, 


l2o( — 
\120 


/ X \ 

\X M ) 


+ 2.1201 — 

x. 


B.A.C.T.— 


= 120 sin 4 1. 


(3.78) 
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For equal coefficients of the normalized variables 
2x m = 120 , i.e. x m = 60. 

2. The computer circuit for solving Eq. (3.78) will consist of the 
usual circuit for a second order equation together with the circuit 
for producing the sin At term. 

3. The scale factors are 


for x, K x = Use ^y* 

for x, K 2 = 6 *o. Use -fa. 

for x, A 3 = i 20 ‘ Use j 2 o’ - 

4. The scaled equation is 

ic _ _ x _ 25/ x \ 
120 ” ~60 ~~ 12\25 ) 


4 * sin At. 


(3.79) 


5. Integrating this gives 



The machine equations are 



and also 

25 25 Jo 60 ^25; 0 

The body starts from rest in the equilibrium position so that 


(3.80) 

(3.81) 

(3.82) 



Fig. 3.14 
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To obtain the desired gains the potentiometer settings are as follows: 

Potentiometer 1 2 3 4 5 6 

Setting 0*2 0*2 0-42 0*24 0*4 0*4 

7. Perform the static check; confirm that the sine wave has the 
correct period. 

8 . Record the values of x/25. The analytical solution is 

x = \6e~ f sin (3* + a).— 12 cos (4 1 — a) 
where tan a = 3/4. 

3.13 Experiment 14. Coupled Circuits 

A steady e.m.f. of 12 volts is applied at time t = 0 to the circuit 
shown. The currents in the two branches after t seconds are denoted 
by x and y and R is 2 ohms. Obtain the graphs of x against t, y 
against t for the ratio of RjL equal to 5, 1, 0*5. 



1. Applying Kirchhoff’s Laws to this circuit the current flowing 
in the arm BC is x — y and the following equations are obtained 


Lx + 2x + 2(x — y) — 12 
Ly + 2y — 2(x — y) = 0. 


(3.83) 


2. The maximum values of the currents x and y may be estimated 
from a consideration of the physical principles involved in this 
circuit. The currents reach a maximum value, i.e. x = 0 = y, and 
hence 


4x — 2y = 12 
4y — 2x = 0 


(3.84) 


giving x = 4, y = 2. 

3. The scale factors are 


for x, 1 /4 
for y, 1/2 
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for x, ~ (see (3.13)) 

for,, 

4. The scaled equations are 



5. Integrating gives 



Initially the currents x and y are zero. Thus 



6 . The computer circuit is 



(3.85) 

(3.86) 

(3.87) 

(3.88) 

(3.89) 

(3.90) 


Fig. 3.16 



ORDINARY DIFFERENTIAL EQUATIONS 41 

7. To obtain desired gains the potentiometer settings are as follows: 

Potentiometer 1 2 3 4 5 

Setting 4/10 L 3/10L 1/10L 4/10L 4/10L 

8 . Perform the static check and then record the values of xj4 
and yj2. 

9. The analytical solution is 

x — 4 — 3e~ 2t/L - e~* t/L 
y = 2 - 3e~ 2tjL + e~ 6t/L 

3.14 Experiment 15. The Force Developed by a Hawser 

A hawser is wrapped around a capstan and a constant force of 200 
lb wt. is exerted at its free end. The force developed at the other end 
depends upon the amount of contact it makes with the capstan and 
upon the coefficient of friction. Examine the relationship between this 
force for various angles of contact up to 4n radians (i.e. two full 
turns) and for coefficients of friction 0-25, 0-2, 0T, 0-02 and 0-01. 



The defining equation of this type of system is found as follows. 
Consider the forces acting in a small element PQ of length 5s. The 
tangential forces due to the tension in the rope are T and T + dT 
and a frictional force F 5s. The normal force acting on PQ is R 5s. 
In the position of limiting equilibrium resolving in the direction of 
the tangent and normal at P respectively the following equations are 
obtained 

(T + ST) cos 59 - T- FSs = 0 ) 

(T + 5T) sin 59 - R 5s = 0 J 

which, as 59 and 5T are small, reduce to 

5T-F5s = 0 | 

T59-R5s = 0 J 

giving F — dTjds , R = T d9jds in the limit as 59 0. 


(3.92) 
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But 


or 


F — jiR 
dT _ dO 

— - jiT — 
ds ds 



(3.93) 


1. This equation has a solution of the form 

T = T 0 e** where T 0 — 200 (see Experiment 7) 


and from this solution the maximum values of T, T m , may be found. 
When 6 = 4n these are approximately 


(a) when ft = 0-25, 

(b) when fi — 0-20, 

(c) when \i = 0T0, 

( d ) when fi = 0*02, 

(e) when jx = 0*01, 


r w = 46 X 10 2 
r m = 25 x 10 2 
r m = 7 x io 2 
r m = 2*6 x 10 2 
r m = 2*3 X 10 2 . 


2. Examination of these values shows that whilst all solutions 
may be recorded using the same scale factor greater accuracy would 
be obtained using a different scale factor in cases (a) and ( b ) from that 
in (c), (( d ) and (e). 

In the former two cases a suitable scale factor would be 10“ 2 /50 
and in the other three 10 _2 /10. 


3. As T 


r 

fxT dO an analogous equation may be set up on the 


Jo 

computer by allowing one second of time in the computer solution 
to correspond to 1 radian of angle in the problem. The equation is 
thus 


T = 


i: 


fiT dt. 


For (a) and ( b ) the scaled equation is 


10 _2 r 

50 



/ lO -2 !’^ 

V"^o 


(3.94) 


For (c), (cl) and (e) the scaled equation is 


o 


(3.95) 
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4. The computer circuit for both scalings is 



Fig. 3.18 


For (a) and ( b ) the output of integrator 1 is 10 _2 r/50, for (c), (d) 
and (e) the output is 10“ 3 T. 

5. To obtain the desired gains the potentiometer settings are as 
follows. 


Case 

Pot. 1 

Pot. 2 

(«) 

0*25 

0-04 

ib) 

0*2 

004 

(c) 

0*1 

0*2 

id) 

0*02 

0*2 

(e) 

0*01 

0*2 


6. Perform the static check and then record the values of T 
against 9. 

7. Check the solutions using T = 2OOe /i0 where (i assumes each of 
the values given in the problem. 

3.15 Further Exercises 

1. Set up a circuit to solve the equation 


with initial condition x = 1 at / = 0, and where a is in the range 
0 ^ a ^ 10. 

[Solution: x = e~ at .] 

2. Set up a circuit to solve the equation 

^I + 9y = 10e-°‘ 
dr 

with initial conditions y = 0, dy\dt = 50 when t = 0. Examine the 
solution for values of a equal to 0, 1 and 10. 
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Solution: v = J— + 

l°a | 

13 

3(a 2 + 9)J 


sin 3 1 
10 


(a 2 + 9) 


cos 3 1 + 


10 


(a 2 + 9) 


-] 


3. Given the simultaneous differential equations 


3x + 2x ~h y = 1 
x + 4y + 3y = 0, 


with the initial conditions x = 0 = y at ? = 0, use the methods of 
simplification of Experiment 12 to show that the greatest expected 
numerical values of x and y are £ and 

Hence program the equations for solution on the computer, set 
up the corresponding circuit and use the static check procedure 
before obtaining the graphs of x and y against ?. 




Solution: 


. * = i - 


y = 


\e — 
— \e 


i 

- 6f/l 1 


-6f/l 


1 ** 


4. Obtain the graphical solution of the equation 


x + 25x = cos 5t 


with initial conditions x = 0-1, x = 0 when t = 0, over an interval 
0 g t <; 10 s. This solution illustrates the case of resonance, the 
amplitude of x increasing with t . The methods of estimating the 
problem variables described previously are not applicable; indeed no 
simple method exists for this type of problem. Assume IxL = 1 and 
|*|* = 5 . 

[Solution: x = 0T (t sin 5 1 -f cos 5/).] 

5. The phenomenon of Beating occurs when two sinusoidal 
signals of similar frequency are added together, resulting in a signal 
of approximately the same frequency but modulated at half the 
difference frequency. Demonstrate this by generating 

x — 0-3 cos 10? and y = 0 3 cos 9 1 
and then producing their sum 


z — 0-3(cos 10? + cos 9?) 

= 0*6 cos 9*5? cos 0*5?. 

Treat x and y as solutions of differential equations as shown in 
Experiment 9. 





CHAPTER 4 


Time Scaling 


4.1 Introduction 

The solutions to the wide range of problems covered by analogue 
computers may extend over periods of time ranging from a few micro¬ 
seconds to many hours. In all cases except those involving a solution 
time of only a few seconds, time scaling of the problem is essential. 

The speed of response of the computer and its associated equip¬ 
ment dictates the minimum solution time taken by the machine. The 
maximum time is essentially determined by convenience in the use 
of the computer and associated recording equipment. 

The need for time scaling is indicated in the scaled equations which 
are obtained as explained in the previous chapter. From these equa¬ 
tions it can be seen that the coefficients of the scaled variables are 
the net gain factors of the inputs to the appropriate integrators. In 
those computers fitted with fixed input resistors and feedback capaci¬ 
tors there is an obvious maximum gain (usually 10) which cannot be 
exceeded, but even where a wide choice of resistors and capacitors 
is available it is acceptable computer practice to work with gains in 
the range 1-50. Hence, if in the scaled equations gains outside the 
acceptable range are indicated time scaling must be applied. 

Consider the following example 

x = 50,000 - 50x - 40,000x (4.1) 

where x = 0 = x initially. 

It can be shown that 0 < x < 10 and -500 < x < 500. 

The scale factors are (i) for x 
(ii) -iV for x. 


The scaled equations are 



(4.2) 

(4.3) 
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Gains of 50, 800, 100 and 50 are required and thus time scaling is 
necessary. 

It can also be seen from this problem that 

a> = 200(see Eq. (3.7)). 

Such a high frequency is likely to produce serious errors in record¬ 
ing equipment and also in other computing components. 


4.2 Methods of Time Scaling 

Time scaling may be performed in two ways. The first of these in¬ 
volves rewriting the problem equation using the following procedure. 
If T = computer time and t = problem time, let 


T — at . 


[Note: If a > 1 the solution is slowed down by a factor a. 
If a < 1 the solution is speeded up by a factor a.] 
The derivatives become 


and, in general, 


d _ d __ d 

dt ~ d{Tja) ~ ^dT 

d 2 d ( d\ 2 d 2 

—= = — — a— I = or —- 

dt 2 d(T/ a)\ dTJ dT 2 


d l = d l 

df a dT"' 


(4.4) 


(4.5) 

(4.6) 

(4.7) 


On slowing down the problem stated in (4.1) by a factor of 100 
the equation now becomes 


10 4 - 5 X 10 4 - 5 x 10 3 — - 4 x 10 4 x, (4.8) 

dT 2 


i.e. 


dT 

x = 5 — 0-5x — 4x 9 


(4.9) 


the dots now denoting differentiation with respect to T. The ranges 
of the problem variables are now 

0<x<10, — 5 < x < 5 

and the new scale factors are 


(i) y for x 

(ii) for x. 
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The scaled equations are 



(4.10) 

(4.11) 


The gains now required are thus 0-5, 8, 1 and 0-5. From Eq. (4.9), 
co = 2^/(63/64) which checks that the problem has been slowed 
by the desired amount. 

The second method of time scaling is carried out by changing the 
rate of integration. This method is preferable to the former as the 
two scaling operations are kept separate and it also reduces time 
scaling to a simple mechanical operation. 

Consider the integration 



■ —11 o 

e { t = \0e h 

_ Jo 


(4.12) 


where e { is a constant. 

If the independent variable t is changed to T by the relation 
T = 5?, then on substitution in Eq. (4.12) 



(4.13) 


These integrations may be represented graphically (Fig. 4.1). 



Comparing Eqs. (4.12) and (4.13) the integrals are seen to be of 
identical form except that the latter is multiplied by a constant of 
In a computer an integrator performs the operation 


e 0 



e t dt + E 0 . 


(4.14) 
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Substituting the new independent variable T this gives 


1 f“' 1 

— e,- dT + E 0 

KCjo a 

(4.15) 

1 f a(1 

e t dT + Eq. 

aiJCjo 

(4.16) 


Thus all that is necessary to change the time scale by a factor of a 
is to change the integrator time constant from RC to olRC . 

As has been shown, solving a differential equation on an analogue 
computer involves successive integrations. The time scaling in any 
such equation may thus be changed by altering the integrator time 
constants as stated. Note also that terms in dxjdt and d 2 xjdt 2 must 
be replaced by a dx/dT and a 2 d z x\dT 2 . 

In interpreting the results the new time scale factor should of 
course be borne in mind. The most convenient way of doing this is 
to relabel the time axis in terms of Tjcc . 

4.3 Summary of Method of Solution of a Differential Equation 

At this stage it is convenient to summarize the various steps involved 
in solving a differential equation. 

(1) Estimate the maximum magnitudes of the problem variables 
and their derivatives. 

(2) Determine the appropriate scale factors and write down the 
scaled equations. 

(3) Draw the computer circuit. 

(4) If time scaling is found to be necessary change the time 
constant (RC) of all integrators to the desired value a RC. 

(5) Apply the check procedures. 

(6) Record the solution and relabel the time axis according to the 
t = Tjoc relationship. 

4.4 Experiment 16. Slowing Down a Second Order 
Differential Equation 

A vibratory system is described by the equation 

x + 2000 x + 10 7 x = 0 (4.17) 

with initial conditions x = 2 and x = 0. Slow down by a factor of 
10 3 and solve. 

1. The system is damped oscillatory and thus the maximum value 
of x may be taken to be approximately 2. Using relations (3.16 and 
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3.17) the maximum values of |x| and |x| are about 6 x 10 3 and 
18 x 10 s . Hence suitable scale factors are 

(a) 1 /2 for x 

(b) 1/10 4 for x 

(c) 1/2 x 10 1 for x. 


2. As 

x = - J* [2000* + 10 7 *] dt 
the scaled equations are 



and 


x 

2 


Jo 


10 4 

~2 


x 

10 4 


dt + 




(4.18) 


(4.19) 

(4.20) 


3. It can be seen from these equations that gains of 2000 and 
5000 would be required. As has been stated, this indicates the need 
for time scaling. Slowing down by a factor of 10 3 the modified 
equations will now be 


x + 2x + lOx = 0 


X 

10 


and 


+ 2 |? 


2 Jo A 10 / 


2 

dT + 


dT + 


X 

To ) 0 




4. The computer diagram is 


(4.21) 

(4.22) 

(4.23) 



Fig. 4.2 
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The potentiometer settings are 

Potentiometer 12 3 

Setting 0*2 0*2 0*5 

5. Perform the static check. From Eq. (4.21) verify that the fre¬ 
quency is 10~ 3 times the frequency in Eq. (4.17). Record the solu¬ 
tion. 


4.5 Experiment 17. Speeding up a First Order Differential 
Equation Given in Problem Form 

Radium decays at a rate proportional to the quantity of radium 
present. In 25 years 0*022 g from a sample of 2 g is found to have 
decomposed. Produce a curve showing the rate of decay to the point 
when one-quarter of the original quantity is left. 

1. If n = number of grammes of radium at any time t years, then 

— = -Kn. (4.24) 

dt 


This gives a solution of the form n = n 0 e~ Kt (see Experiment 7). 
It is given that n 0 = 2 and when t = 25, n = 1-978. Hence 
1-978 = 2e~ K - 25 , i.e. log c 0-989 = -25K giving K = 0-00044. 

2. The maximum value of n is 2 and an exactly similar computer 
equation is obtained if 1 s of time in the computer is made to cor¬ 
respond to 1 year in the problem. The describing equation is thus 

— = —0-00044«. (4.25) 

dt 


The scale factor for n is 

3. As 


n = 

the scaled equation is 


m t 

0*00044 n dt + ( n) 0 

Jo 


n 

2 


0*00044 




dt + 



(4.26) 


(4.27) 


Speeding up by a factor of 10 3 (i.e. 1 s on computer = 10 3 years) 
a modified equation is obtained: 



(4.28) 
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4. The computer circuit for this is 



The potentiometer settings are 

Potentiometer 1 2 

Setting 0-44 1 

For a better recording a different time scale factor may be necessary. 
Repeat the experiment speeding up Eq. (4.25) by factors of 100, 
10 , 000 . 


4.6 Experiment 18. Xenon Poisoning in a Nuclear Reactor 

A reactor which has been operating at a high thermal-neutron flux 
level is suddenly shut down. Iodine 135 and xenon 135 are then 
present in amounts 3-17 x 10 15 and 1-90 x 10 15 atoms per unit 
volume. The half-life periods of the two nuclides are 24,120 and 
33,120 hr respectively and it is known that iodine 135 decays to 
xenon 135 which in turn decays to comparatively stable caesium 135. 
Obtain graphs showing the concentrations of iodine 135 and xenon 
135 as functions of time after shut down. 

1. It can be shown that the equation for the concentration of 
xenon is 

^ = -X t X + X 2 I 
at 

and that for the concentration of iodine 


where and a 2 are the decay constants per second for xenon and 
iodine respectively, and X, I are the concentrations of xenon and 
iodine respectively in atoms per unit volume. The solution of the 
first order equation dy/dt = —Ky is y = y 0 e~ Kt , and applying I his 
to the problem, the describing equations are found to be 


dl 

dt 


0-6931 


I = 


-2-874 x 10“ 5 / 


24,120 


(4 29) 
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— = _0^931 x + = _ 2 . 092 x 10 -s z + 2 . 874 x 10 -sj 

dt 33,120 24,120 

(4.30) 

2. The maximum value of/is obviously 3-17 x 10 15 . An estimate 
of the maximum value of X may be obtained in a similar manner to 
that used in Experiment 12. It is easier however to try say a maxi¬ 
mum of 4 x 10 15 ; if this is too small overloading will occur but it 
will not be too difficult to re-scale. 

3. Using the above values the scale factors are 

(i) for /, i x 10“ 15 

(ii) for X, i x 10“ 15 . 


4. The machine equations are 


!22!/=-f 2 . 

4 Jo 

l 

-15 [‘I ~ 

—X = - I 209: 


874 x 10 


-YKT 15 


I \ dt 4- 


092 x 10“ 5 i 


-2-874 x 10" 5 


/10~ 15 Y] 


Speeding up by a factor of 3600, i.e. 1 s of computer time equivalent 
to 1 hr of the problem time, these equations reduce to 


10 “ 15 r r /to -15 \ 

= -J 0T03L —^—jj dT + 0-793 (4.33) 



Fig. 4.4 
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6. The potentiometer settings are 

Potentiometer 1 2 3 4 5 

Setting 0 103 0-793 0-075 0-103 0-475 

7. Carry out the check procedure and obtain the graphs of I 
against t and X against t. 

4.7 Experiment 19. Simultaneous Differential Equations 
Requiring Time Scaling 

Obtain the graphical solutions of the equations 

x -f lOOx + 10 4 x + 10 4 y = 4 x 10 4 
y + I00y + 9 x 10 4 y + 200x - 19 x 10 4 x = 22 x 10 4 

given that x = y = 0; x = 200; y = 400 at t = 0. 

1. As in Experiment 12 nothing is known of the maximum values 
of the problem variables, but estimates may be obtained using a 
procedure similar to that outlined in that experiment. 

Consider the second equation with the x and x terms zero. The 
modified equation is 

y + 100y |9x 10 4 y = 22 x 10 4 . (4.36) 

Using the equal coefficient rule the following estimates of maximum 
numerical values are obtained: 

\y\ m = 22 x 10 4 , \y\ m = 22 x 10 2 , \y\ m = 

{Note: These values are calculated on the basis of zero initial condi¬ 
tions and will be only approximate.) 

In the problem use \y\ m = 25 x 10 2 and \y\ m = 10. 

Substituting this maximum value of y in the first equation of 
(4.35) gives 

x + 100 x + 10 4 x =-6 x 10 4 . (4.37) 

From this the estimated maximum values are found to be 
|jc| m = 6 x 10 4 , |x| m = 600, [x| m = 12 by further use of the equal 
coefficient rule. 

Use in the problem |x| m = 10 3 and |x| m = 10. 

2. The machine equations may now be written 



B.A.C.T,— 5 
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3. Slowing down by a factor of 100 the modified equations are 



The dots now denote differentiation with respect to T, where 
T = 100 L 

4. The computer diagram is 



Fig. 4.5 
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5. The potentiometer settings are 


Potentiometer 

i 

2 

3 

4 

Setting 

04 

0-8 

0-76 

0-88 

Potentiometer 

5 

6 

7 

8 

Setting 

0-36 

0-25 

0-2 

016 


6. Carry out the check procedure and then obtain the graphs of 
x against t and y against /. 

7. The theoretical solution is 

x = sin 200 1 — i cos 200/ + i 
y = 2 sin 200/ — 3-5 cos 200/ + 3*5. 

4.8 Experiment 20. A Suspension Problem 

A two-wheeled trailer chassis weighs 960 lbf and its axle 128 lbf. 
The spring constant between the chassis and axle is 800 lbf/ft and 
that between the axle and ground (due to the tyres) is 1500 lbf/ft. 
Examine the effect of various shock-absorber damping coefficients 
upon the motion of the chassis when the trailer mounts curbs of 2, 
4 and 6 in. 

A simplified diagram illustrating the sus¬ 
pension system of one wheel is shown in 
Fig. 4.6. 


M t = lb = 480 lbf 
M 2 = lb = 64 lbf 
K t = 800 lbf/ft 
K 2 = 1500 lbf/ft 

c = shock absorber damping coefficient 
variable from 10 to 100 lbf s/ft 
x 3 = 0T66; 0*333; 0-500 ft. 

1. Comparing this problem with that in Experiment 13 and using 
Newton’s second law the equations of motion for the two masses are 
found to be 

Mi 

— 1 -x 1 + c(x x — x 2 ) + K l (x l — x 2 ) = 0 (4.46) 

9 

—x 2 + c(x 2 - xO 4- K x (x 2 - x x ) + K 2 (x 2 - x 3 ) = 0, (4.47) 
9 


Mx 








m 2 






Fig. 4.6 
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i.e. 

15*! 4- c(x x — x 2 ) 4- 800^! — x 2 ) = 0 (4.48) 

2x 2 4- c(x 2 - Xi) + 800(x 2 - x x ) 4- 1500(x 2 - x 3 ) = 0. (4.49) 
Initially x x = 0 = x 2 and x x — 0 = x 2 . 

2. Practical considerations enable a good estimate of the maxi¬ 
mum values of x x and x 2 to be made in this problem. It will be 
assumed that the maximum displacements will not exceed twice the 
value of the input function x 3 . 

3. Consider Eq. (4.48) without the presence of the terms contain¬ 
ing x 2 and x 2 

15x x 4- cx x + 800x x = 0. (4.50) 

In the case where damping is least c = 10, hence 

15Xi + 10x t 4- 800x! = 0 (4.51) 

and (D n = = 

Then as |x t | m = 1, \x t \ m = IO^/it- The value |*i|m = 5 ft /s will be 
used. 

Substituting these maximum values in Eq. (4.49) 

2x 2 + 10x 2 - 50 + 800x 2 - 800 + 1500x 2 - 1500x 3 = 0, (4.52) 


i.e. 

2x 2 + 10x 2 + 2300x 2 = 850 + 1500x 3 (4.53) 

co n = 107\ 3 - and |x 2 |„ = 1, |* 2 |« = lO^. The value 

|x 2 | m = 50 will be used. 

4. The scale factors are 


for Xi and x 2 , 1 

for x u i 
for x 2 , -5V. 


5. The scaled equations are 


■&)- 


2Jx 2 


3150 


—. Xl - i°°.x 2 "| dt (4.54) 
75 75 


*2 

50 



c( x 2 
2l 50 


c_lx_ 
20V 5 


\ I + 23x 2 - 8xj - 15 x 3 dt (4.55) 




(4.56) 
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* 2= “1o- 50 &H'- (4 - 57) 

All initial conditions are zero since the system starts from rest in the 
equilibrium position. 

These equations may be slowed by a factor of 10 to give 


.1. 


+ ~ “;<ij <iT (4.58) 


*1 _ _ f T T cJxA _ c^(x2 

0-5 Jo |_150\0-5y 15^5 

— = - f + 2-3x 2 - 0-8xj - 1-5 x 3 1 dT 

5 JoL 20 \ 5 / 200y0-5y 1 3 J 


Xl = - f r -H ^ dT 
Jo 2\p-5/ 

-<th 


(4.59) 

(4.60) 

(4.61) 


6. The computer circuit is 



Fig. 4.7 
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7. The potentiometer settings are 
Potentiometer 

1234567 8 9 10 11 

Setting 

c 8 c _ i -L 1-5*3 A. 0-23 0-5 

150 75 150 75 2 200 200 10 

The initial condition settings are all zero. 

8. Carry out the check procedure and for the different values of 
x 3 examine the effect on the chassis vibration of an increasing damp¬ 
ing coefficient c. 

4.9 Further Exercises 

1. Set up a circuit to solve the equation 

x + 100 x + 10 6 x = 10 6 
given that x = 0 = xatt = 0. 

[Approximate solution is x = 1 — e~ 50t sin (10 3 t 4* T52).] 

2. The velocity of a chemical reaction is proportional to the con¬ 
centration of the reacting substance. Show that if a be the initial 
concentration of the reagent and x the amount transformed after 
time t then 

dx 

— = K(> a - x). 
dt 

Set up a suitable circuit to obtain a graph showing the variation of 
x with t over the first second of the reaction when K — 5-2 and 
a = 90-3. 

[Solution: x = 90-3 (1 — e 5 ' 2 /)-] 

3. Examine the solutions of the equation 

x + 0 06 x + 0-011 x = 0-011 sin cot 

for co = (i) 0-08, (ii) 0-09, (iii) 0-1, given that x = 0 = x at / = 0. 
[Solution: 

* = e -°' 03, [A cos0-1005t + B sin 0-1005/] + K sin cot + Lcos cot 
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where, for 


(i) A = 1-195, 

B = -0-5547, 

K= 1-145, 

L = 

-1-195 

(ii) A = 1-581, 

B = -0-2885, 

K = 0-849, 

L = 

-1-581 

(iii) A = 1-783, 

B = +0-2364, 

K = 0-297, 

L = 

-1-783.] 


4. Obtain the graphical solution of the simultaneous equations 

x + 10” 2 x + 10” 4 (x - y) = 0*5 x 10” 4 

y + 10 ~ 2 y + 10“ 4 (x + y) = 1-5 x 10“ 4 

given that x = -l,y = 1*5, x — 10” 2 , y = 2 x 10” 2 at t = 0. 

[ Solution: x = sin 10 ~ 2 t — 2cos 10“ 2 t + 1 "1 
y — 2 sin 10 ~ 2 t + cos 10 ~ 2 t + 0*5.J 






CHAPTER 5 


The Use of Non-Linear Units 


5.1 Introduction 

In the preceding chapters the problems were restricted to those 
involving linear operations only. Terms such as y dy\dx or yjy did 
not occur in the equations to be solved. This limitation may be lifted 
if certain other non-linear elements are included in the computing 
equipment. The analogue computer handles non-linear differential 
equations just as easily as linear ones. 

A brief outline of the most important non-linear unit—the multi¬ 
plier—will be given in this chapter together with some of its major 
applications. Experiments utilizing these techniques will then be 
given. For a thorough review of multiplying units the student should 
consult the computer reference books and the operator’s manual for 
the particular computer at his disposal. 

5.2 The Multiplier 

A computing potentiometer is the simplest form of multiplier. It will 
multiply a varying voltage by a constant less than or equal to unity. 

To form the product of two varying voltages special units have to 
be designed. Many forms of these are available but broadly they may 
be classified as the servo-multiplier type or the electronic multiplier 
type. 

+ 1M.U. +Y 



Servo-multipliers are based on the feedback principle (see Fig. 
5.1). The multiplier consists of a servo-driven follow-up potentio- 
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meter, F\ and several multiplier potentiometers (three, A, B, C, have 
been shown in Fig. 5.1). The wipers of these potentiometers are 
mechanically ganged together, indicated by the dotted line in the 
figure. 

If an input signal X M.U. causes the feedback potentiometer 
wiper to rotate 6 radians, then 


- = X- (5.1) 

n 

Thus if signals + Y and — Y M.U. are applied as shown, the output 
of the wiper of potentiometer A is 

ny 

e 0 = = XY M.U. (5.2) 

n 

The potentiometers B and C may also be supplied with inputs, say 
+ U and — U and + V and — V, to provide further products + UX 
and + VX respectively. Sign inversions may be obtained by reversing 
the polarity of the connections of the appropriate multiplying 
potentiometers. 

The static accuracy of the servo-multiplier is largely governed by 
manufacturing tolerances on the potentiometers, these being sig¬ 
nificantly small for most applications, but the most severe limitation 
is the speed of response. 

The range of frequencies at which the servo-multiplier will operate 
accurately is rather low, although if only one of the functions X , Y 
is of high frequency, this signal may be connected to the function 
potentiometer A and the slower signal used as input to the servo. 

The resistor R is a load compensating resistor, equal in value to the 
load on the function potentiometer. If R were not present the loading 
on the potentiometer A caused by the input resistor of the next 
computing element would cause a small loss in the output voltage 
from the multiplier. By loading the feedback voltage in an identical 
way this type of error can be avoided. 

Other errors may be due to short-comings in the electrical and 
mechanical zeroing of the potentiometers F, A, B , C. 

Four quadrant multiplication is possible, though if it is not 
necessary a sign reversing amplifier can be saved by arranging that 
the servo is driven to the half of the function potentiometer that is 
supplied with the Y signal. 

Electronic multipliers have been designed in many ingenious 
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forms, some of the more familiar types being the Quarter Square, the 
Time-Division, the Logarithmic, the Photomultiplier, and the Hall 
Effect. Basically each will give an output e 0 = — e x e 2 M.U. where 
e x and e 2 are the inputs. Usually the output is developed from a 
computing amplifier having very low output impedance, while the 
input impedance varies from one type of multiplier to the next. 

The frequency response is much better than in the case of the 
servo-multiplier and while accuracy was at one time a limiting 
factor modem electronic techniques are now producing some very 
accurate multipliers. 

5.3 The Dividing Circuit 

Division may be performed by inserting a multiplier in the feedback 
loop of an operational amplifier. Figure 5.2 shows how this is done. 



The signal X is fed to an input of a high gain amplifier through a 
resistor R and another signal Y fed to one input of the multiplier. 
If the output of the amplifier, e 0 , is connected to the second multiplier 
input the output will be — e 0 Y from Eq. (5.1). This provides the 
second input to the amplifier through a resistor R. 

Assuming the grid of the high gain amplifier to be at earth 
potential, ji being very large, KirchhofTs first law gives 

— - = 0, (5.3) 

RR 

i.e. 

*o = f • (5.4) 

In this way division of X and Y is achieved. From the last result it 
is evident that Y should be greater than X f or else voltages greater 
than one machine unit would result. However, this can be obviated 
by compensating in the input resistors as in Fig. 5.3. 
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R 2 X 
e 0 ~ —- —. 
R t Y 


(5.5) 


If X exceeds Y then make R t greater than R 2 by the same ratio. 


5.4 The Squaring and Square Rooting Circuits 

Other applications in which the multiplier may be used are those of 
squaring and square rooting. For squaring, the multiplier is used 
with equal inputs X the output being — X 2 (see Fig. 5.4). 



Fig. 5.4 


In square rooting the multiplier is used as a squaring device this 
being placed in the feedback path of an operational amplifier as 
shown in Fig. 5.5. 


X o AWW- 

K 



Fig. 5.5 


Using a similar argument to that for the division circuit 


— — — = 0, 

(5.6) 

R R 

=v*- 

(5.7) 


Squaring devices may be produced other than by using a multiplier. 
For instance a system of biased diodes may be arranged so as to 
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produce a straight line approximation to a square law. This forms the 
basis of most squaring units. 

The generation of a quotient and of a square root are examples 
of the implicit function technique. An implicit equation between 


the input functions and the output function is 
of division the equation solved is 

solved. In the case 

O 

II 

O 

I 

(5.8) 

and for square-rooting 


X - el = 0. 

(5.9) 

More generally the equation 


F(e o, Y) + JT = 0 

(5.10) 


where F is an arbitrary function covers a wide variety of implicit 
functions. Figure 5.6 shows the basic circuit for this equation. 



Other forms of non-linear units are available, such as sine and 
cosine resolvers, which perform these operations on any input volt¬ 
age; diode limiting units; backlash and hysteresis units; relay 
comparators; and arbitrary function generators. Although experi¬ 
ments involving the use of these units are beyond the intended scope 
of this book they nevertheless are extensively employed in many 
branches of analogue computing. 

5.5 Experiment 21. The Use of a Multiplier and Divider 

The binomial expansion of 1/(1 + x) is 1 — x + x 2 — x 3 + 
Obtain graphs of 1/(1 + x ) and the first four terms in its expansion, 
and hence show that errors of more than about 2 per cent of full 
scale occur when these four terms are used to represent 1/(1 + x) 
when x is greater than about 0-35. 

1. The function x can be obtained as a linearly increasing signal 
by an integrator with a constant input. 

To allow a reasonable computing time for recording purposes, sup¬ 
pose that a 10 s run is required. Then x must reach 1 M.U. after 10 s. 
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Thus let 

. __ t 

x ~io 

where t is problem time. 

Two multipliers are used to provide the signals x 2 and x 3 . Further 
sign reversing amplifiers may be needed depending upon the type of 
multiplier used. 

2. The basic computer diagram is shown in Fig. 5.7. 


-1 M.U. 



3. Record the output of amplifier 2 for a period of 10 s full-scale 
calibration. Check the curve for accuracy using suitable values of x 
in the range 0 x ^ 1. 

4. The function 1/(1 + x) can be generated using the method of 
division already discussed, where the divisor is 1 -f x and the divi¬ 
dend constant equal to unity, 

5. The basic computer circuit is shown in Fig. 5.8. 



Potentiometer 12 3 

Setting 005 0-5 0*5 

In order that the output of amplifier 1 does not exceed 1 M.U. a 
signal ^(1 + x) is generated. To compensate for this the numerator 
is made 0*5 using potentiometer 3. 



66 


BASIC ANALOGUE COMPUTER TECHNIQUES 


6. Record the output from amplifier 2 for a period of 10 s on the 
same scale as the previous graph. Check the curve by calculation 
using suitable values of x. 

7. Compare the two curves and determine the point at which the 
difference is 2 per cent of full scale. 


5.6 Experiment 22. A Further Example on Division 

The purpose of this example is to outline a method which may be 
used when the divisor is zero at time t — 0. 

Consider the function y = (i sin t)jt. Obviously when t = 0 
overloading might occur in the division circuit. This difficulty is 
minimized by generating a function of the form Ke~ at and adding 
this to the divisor. When t — 0 the divisor is of magnitude A and the 
exponential function diminishes at a rate depending upon the value 
of a. It is usual to adjust a so that a voltage is provided until the time 
divisor is built up. (Note: in a time of 4/a seconds the exponential 
output will fall to under 2 per cent of its initial value.) The value of K 
is made just large enough to prevent 
overloading. 

1. In order to obtain a reasonable 
length of computing time generate ^ 
sin t and t\ 10. 

2. The extra term in the divisor (in 
this case O-Ole” 10f ) is formed from the 
circuit of Fig. 5.9. 

3. The complete circuit is 
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Potentiometer 12 3 

Setting 005 0-1 0*01 

5.7 Experiment 23. Use of a Multiplier to Perform the 
Square-Root Operation 

When a heavy cable (IFlbf per unit length) is hanging in equilibrium 
under the action of a horizontal force H lbf as indicated in Fig. 5.11 
it can be shown that 




Assuming that the maximum value of y is 20 ft and that of dy\dx 
is 2 obtain the graphs of y against x for values of 


W 

H 


i . 
1 0 0 > 


1 . 
2 0 0 ’ 


1 

3 00 * 


The values of y and dyjdx when x = 0 are both zero. 

1. If x = t then the describing equation may be written 

y = + y 2 )- 

H 

2. The scale factors are 

for y, 2 ^ 
for y, 

3. The scaled equations are 



where a = WjH. 


(5.12) 


(5.13) 

(5.14) 
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Speeding up by a factor of 10 these now become 



where the dots denote differentiation with respect to T. 
4. The computer diagram is 


(5.15) 

(5.16) 



Fig. 5.12 


Potentiometer 1 2 

Setting 10a 0-25 

5. The theoretical solution of Eq. (5.11) is 


H( U W 1 . 
— cosh —x — 1 |. 
W\ H 


5.8 Experiment 24. Solution of Mathieus’ Equation 

This equation describes the behaviour of many physical systems 
including wave guides, frequency modulation and sinusoidally 
excited mechanical systems. The equation may be expressed in many 
different forms one of which is 

^4- + (a — 2b cos wt)y = 0. (5.17) 

dt 2 

The case in which a = 2b, w = 2 and 0 ^ a ^ 5 provides some 
interesting studies of stability and instability. In the equation the 
initial conditions are taken to be y = 1 and y = 0. 
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1. The equation now becomes 

y + a{\ — cos2f)j — 0 (5.18) 

from which it can be seen that it is necessary to generate a function 
x = 1 — 2 cos 2 t. An accurate method of doing this is by solving a 
differential equation, thus if 


x = \ — cos 2 1 

(5.19) 

then 


x = 2 sin 2 1 

(5.20) 

and 


x = 4 cos It = 4(1 — x). 

(5.21) 

As the maximum numerical values of x and x are 2 the scale factors 


are ^ in both cases. 

The scaled equations are 


and 


2 JoV 2 

X = _f‘ _X 

2 Jo 2 


dt + l- 


dt + 


(5.22) 


(5.23) 


2. As the maximum value of x is 2 the maximum values of y and 
y may be estimated from an examination of the equation 


y + lay = 0, 


(5.24) 


y = 1, y = 0 at t = 0. This has a solution of the form 

y — cos yjlat (5.25) 

whence the maximum value of y is 1 and that of y is yj 10. However, 
as has been stated, some of the solutions of the problem are unstable. 
A larger maximum value of y is therefore advisable. Let this maxi¬ 
mum value be 5. The maximum value of y will now be 5^/10. Take 
this as 20. The scale factors for y and y are | and respectively. 

3. The scaled equations are 



(5.26) 

(5.27) 


B.A.C.T.— 6 
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4. The computer diagram is 



5. The potentiometer settings are 

Potentiometer: 1 2 3 4 5 

Setting: — 0-4 0-2 0-4 0-2 

20 

6. Perform the static check and obtain the graphical results 
which are as shown in Fig. 5.14. 


5.9 Further Exercises 

1. Water flows through a circular orifice of radius r ft from the 
side of a cylindrical tank of radius R ft. The velocity of flow is given 
by v = *J(2gh)ftl$ where h ft is the height of the surface of the water 
in the tank above the centre of the orifice. 

Show that the equation describing the rate at which the surface 
is falling is 

l - 

If R = 5 and h = 100 initially, show how h varies with t for 
r = 01, 0-2, 0-3, 0 4. 

|^The theoretical solution is yj'h = 10 — , taking g = 32.J 

2. Set up a suitable circuit to examine the variation of v with t 
for the equation dvjdt + Av 3 — w 2 r = 0 where 


(i) 


"I/ 


dt. 


(ii) w = 2 x 10 4 , 

(iii) 0 < r < 1, 


(iv) at t = 0, v = 0 and r = 0-5, 

(v) Max. value of v = 2 x 10 4 , 

(vi) 0-75 x 10~ 4 < A < 2-5 x 10“ 4 . 
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By varying A within the defined limits obtain the type of v-t curves 
illustrated in Fig. 5.15. 



Fig. 5.15 


3. An important equation occurring in the theory of feedback 
oscillators has the form 


y + viy 2 - l)y + y = 0. 


This is known as Van der Pol’s equation. For small (but arbitrary) 
initial values of y and y the steady state solution is a periodic wave 
form, varying in character with different values of jjl. Assuming the 
maximum numerical values of y and y to be 2 and 4 respectively, 
show that the machine equations are 



Obtain graphs of y against t and y against y [an X-Y recorder is 
essential in this case] for values of n equal to 0, 0*4, 2. In these 
experiments use (y) 0 = 0 and (y) 0 = 0 02. In the case of the phase 
plane plot (i.e. y against y) note that the solution is always limited 
by a definite bounding curve. This is known as the Limit Cycle and 
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is of great importance in the examination of the stability of oscillatory 
circuits. For a given value of \i the same limit cycle is obtained ir¬ 
respective of the initial condition values. Typical results are as shown 
in Fig. 5.16. 




CHAPTER 6 


The Use of Transfer Functions 

6.1 Transfer Functions 

In the study of equations representing physical systems, it is often 
convenient to consider the relations among the input and output 
variables of the system. For linear systems these relations are ex¬ 
pressed as functions known as Transfer Functions. 

Consider, for instance, a stable system governed by the equation 

a ~r£ + b i 7 + c * =/(0 (6.i) 

at at 

where bja is positive, and assume all initial conditions to be zero. 

Taking Laplace transforms of both sides of Eq. (6.1) produces the 
algebraic equation 

(as 2 + bs + c)x(s) = f(s) (6.2) 

where s is the Laplace transform variable, in general complex, and 
x(y) and f(s) are the Laplace transforms of x(t ) and f(t) respectively. 

Students unfamiliar with the technique of the Laplace transforma¬ 
tion are referred to any standard mathematical text on this subject. 
An adequate account may be found in Chapter 1 of Analog Computa¬ 
tion by A. S. Jackson. 

The function 

as 2 + bs + c 

is called the characteristic function , and the equation 

as 2 + bs + c = 0 (6.3) 

the characteristic equation (corresponding to the homogeneous case 
of Eq. (6.1)). The roots of Eq. (6.3), which may be complex, deter¬ 
mine the form of the complementary functions in the solution of 
Eq. (6.1). For example, if the two roots are s — s t and s = s 2 the 
complementary function is of the form Ae Slt + Be S2 \ A and B being 
constants. 
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Since Eq. (6.3) is algebraic in s it yields the ratio 

^ = —-— -. (6.4) 

/(s) as 2 + bs + c 

The function on the right-hand side is the transfer function relating 
x(s) and f(s) and is the ratio of the Laplace transform of the output 
of the system divided by the Laplace transform of the input (or 
driving function). Throughout this chapter, Y(s ) will be used to 
denote a transfer function. In the case of Eq. (6.4), 

y(s) = (6.5) 

m 

Thus Eq. (6.1) may be represented as shown in Fig. 6.1 by a block 
corresponding to the transfer function, with f(s) as input and x(s) 
as output. 


f(s) 

n , 1 

x(s) 


rw as 2 + bs + c 



Fig. 6.1 


It is to be noted that a transfer function gives no account of the 
initial conditions of an equation; these were assumed to be zero 
in the formulation of the transfer function. 

When the driving function f(t) is sinusoidal, the particular integral 
of Eq. (6.1) is also found to be sinusoidal, having the same frequency 
as f{t) but differing in phase and amplitude from it. If, for instance, 

fit) = C sin cot (6.6) 

then the particular integral takes the form 

x(t) = D sin (cot -f- 0) (6.7) 

where (j) is the phase difference between f(t) and x(t) and C and D 
are constants. Other forms of input may usually be expressed as a 
Fourier series of sine and cosine terms. Hence it is important to study 
the behaviour of a system when it is subjected to sinusoidal inputs. 

In the solution of an equation representing a stable system the 
particular integral is essentially the only part left after a large time 
t has elapsed; for this reason it is termed the steady state solution . 

One way of measuring the performance of a system is to examine 
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the steady state solution of the describing equation for driving 
sinusoids covering a whole range of frequencies. In other words the 
dependence of D and 0 on co is observed. This is known as the 
sinusoidal steady state frequency response, or more usually, the 
frequency response. 

Mathematical convenience is achieved by using a complex driving 
function of the form 

fit) = Ce jcot (6.8) 

= C(cos cot + j sin cot), (6.9) 

where j — — 1. 

The imaginary part of this function is C sin cot and the resulting 
complex steady state solution is He jo3t y where 77 is complex and 

He* m = |ff|{cos (cot -f arg 77) + j sin (cot + arg 77)}. (6.10) 

Comparison of the imaginary part of this solution with the right- 
hand side of Eq. (6.7) gives 

D-\H\ (6.11) 

and 

(j) = arg 77. (6.12) 

Consequently, the variation of D and 0 with co can be established 
from the complex steady state solution. 

If the Laplace transform solution of Eq. (6.1) is pursued with 
f(t ) = Ce i( ° t 9 the resulting steady state solution, x ss (t), is found to be 


xjt) = 

Ce Ja> ‘ 

(6.13) 

a (/to) 2 + bjco + c 

xjt) _ 

1 

(6.14) 

Ce Ja “ 

a(jco) 2 + bjco + c 


Equation (6.14) determines how the ratio 

steady state solution 
driving function 


varies with frequency. Clearly, this may be obtained by replacing 5 
by jco in the transfer function Y(s), giving 


Y(jco) = 


_ 1 _ 

a(jco) 2 + bjco + c 


(6.15) 
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This is called the frequency-domain transfer function . This may be 
written in the exponential form 


70) = G(co)e j<f,M 

where 

G(<w) = |YO)( 

_ magnitude of sinusoidal output 
magnitude of sinusoidal input 
and 


(6.16) 


(6.17) 


<0) = ar E Y(ja>) 

= (phase angle of output) — (phase angle of input). (6.18) 

G(co) and (j)(co ) are the gain and phase characteristics respectively 
of the transfer function. Consideration of these at various angular 
frequencies a> will determine the frequency response. 

Although the example quoted above deals with a second order 
differential equation, transfer functions are obtained in the same 
way for linear differential equations of any order and the definitions 
of G(co) and (j)(a)) given in Eqs. (6.17) and (6.18) remain valid. In 
each case the transfer function 7 ( 5 ) is the ratio x(s)ff(s). 


6.2 Simple Block Diagrams 

As mentioned earlier, a system governed by a differential equation 
may be represented diagrammatically as in Fig. 6.1 by a block, or 
black box, consisting of the transfer function Y(s) with input J(s) 
and output *($), the transforms of the input (driving) function and 
the output function (solution) respectively. More complex systems 
are usually divided into several distinct sections, individually de¬ 
scribed by a differential equation of reasonably low order, and each 
represented by a transfer function. The product of the individual 
transfer functions is the overall transfer function. For instance, if two 
transfer functions, Yfs) and 7 2 (^), are linked to their corresponding 
input and output transforms by the equations 

(6.19) 

and 

x 2 (s) = r 2 (s).x 1 (s) (6.20) 

then the resulting equation relating the input transform f t (s) and 
the output transform x 2 (s) is 

x 2 (s) = Y 1 (s)Y 2 (s)f l (s). 


( 6 . 21 ) 
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Hence 

= n*) = n(s)n(s). (6.22) 

JM 

Any number of transfer functions may be multiplied in this way. 
Schematically the black boxes representing the various transfer 
functions are connected in cascade as in Fig. 6.2. For linear systems, 



Fig. 6.2 


the order in which Y^s) and y 2 (.?) occur is not important so long as 
the signal represented by x x (s) is not required. In connecting blocks 
in this way it is important that one block draws no energy from the 
preceding one. As far as actual computer circuits are concerned, this 
implies that electrical loading does not occur. In practice the loading 
is usually kept down to negligible proportions. 

In systems which employ feedback* the output may be passed 
through a further block to an error sensing device—usually the 
input grid of an operational amplifier and compared with the input 
function. Figure 6.3 shows the block diagram for a simple feedback 
system. 



Fig. 6.3 


Y f (s), Y b (s ) are the transfer functions of the forward path and 
feedback path respectively, and the error device C algebraically 
sums the inputs to it with the signs shown; that is 



£ = /(*) - x 2 (s). 

(6.23) 

But 


X 2 (s) = Y^S^X^s) 

(6.24) 

and 


a x 0) = Y f (s)e. 

(6.25) 
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This gives the overall output/input relation 
xi(s) _ Y f (s) 


M 1 + Y f (s)Y b (s) 


(6.26) 


Y f (s ), Y r (s)l[ 1 + Y f (s) F 6 (.y)] are termed the open loop and closed loop 
transfer functions respectively between x ± (s) and f(s ). 


6,3 The Use of RC Networks to Simulate Transfer Functions 

Many methods are available for setting up transfer functions on a 
computer. One of these is the use of complex RC networks as the 
input and feedback impedances of an operational amplifier. 

Let Z/,Zi represent the complex impedances in the feedback and 
input paths of Fig. 6.4. As before, assuming an earth potential 



exists at the amplifier grid, it may be shown that 


foOO = _ Z f (s) 
e^s) Z,(s) 


(6.27) 


where e 0 (s) and e t (s) are the Laplace transforms of the output and 
input voltages respectively. 

Thus 

_ § 7 ? sF(5) - (6 - 28) 

Z&s) 


The impedances of a resistance R and capacitance C are R and 
l/Cs respectively. Hence, for the circuit of Fig. 6.5 


R 

I—VWW-1 
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or 


y(s) = 


R 

ijCs 


= -RCs 


e 0 (s) = - RCse^s). 


On taking inverse transforms 

e 0 (t) = -RC 

at 


(6.29) 

(6.30) 

(6.31) 


The transfer function is — RCs and the operation performed is 
pure differentiation with gain ~RC. This result may be arrived at 
using the elementary method given in Chapter 1 for the integrator. 

The differentiating circuit, as it is shown in Fig. 6.5, is not of very 
great practical use, since, if an unwanted noise signal of high fre¬ 
quency content is present along with the true signal, the output of the 
differentiator inherits this noise multiplied in magnitude by its 
angular frequency. The resulting pen recorder trace may exhibit a 
large amplitude jitter superimposed on the correct solution. In addi¬ 
tion, amplifiers connected in this way are liable to run into satura¬ 
tion and instability. 

Methods of minimizing these difficulties are based on the filtering 
off of the high frequency components of the signal and lead to 
approximate forms of differentiation. Details of this can be found in 
Experiment 25. 

Another commonly occurring circuit is that of Fig. 6.6. This 


Ri 



represents the simple lag transfer function 


7 ( 5 ) = 


G 

1 + Ts 


(6.32) 


where G = R 2 /R i and T = R 2 C. 
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The result is easily obtained using Eq. (6.27) 


Z /= l 


z, = Ri 


— + Cs ) = - — - 

R-2 J 1 + R. 2 CS 


Y(s) = - 


^ 2/^1 
1 -|- R 2 GS 


(6.33) 


T is called the lag time constant in seconds, and the transfer func¬ 
tion represents the output/input relation between two variables 
e 0 , satisfying the first order differential equation 


T 


de 0 

dt 


+ e o 


-Ge t . 


(6.34) 


Other simple transfer functions and their circuits are shown in 
Fig. 6.7. 


6.4 Checking the Accuracy of Transfer Functions 

The use of transfer functions in analogue computation involves just 
as much caieful checking as do the normal methods for solving 
differential equations. These checks may be carried out by ascertain¬ 
ing that the gain and phase curves for each transfer function block 
are correct. 

Measurements are taken of G(co) and (f>(co) for sinusoidal inputs 
of known character over a wide range of frequency, and the results 
compared graphically with the calculated values. The two should 
agree to within the degree of accuracy involved in the computing and 
measuring equipment. 

It is customary to plot G in decibels (i.e. 20 log 10 G) and <£ in 
degrees on linear scales against angular frequency co on a logarithmic 
scale. The plots are termed Bode Plots. 

Another method of presenting the same information is to plot 
G against <£ in polar form—called a Nyquist Plot—in which case 
margins of stability in gain and phase are more easily abstracted. 

Other methods of testing transfer functions include the transient 
response and impulse response techniques. 

In the former a step function is applied to the input of the transfer 
function and the output/time graph recorded and compared against 
the calculated, while in the latter the input is an impulse function. 
Essentially these two techniques are used to analyse the system 
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X2 - R2C1 

T3 == R3C2 


Fig. 6.7 
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under test in the time domain, while the Bode and Nyquist plots 
are contained in the complex frequency domain. A correspondence 
may be set up between the frequency response and the impulse 
response and vice versa so that information on the system in one 
domain can be transformed into information in the other. Discussion 
of this is beyond the scope of the present volume. 

The ideas put forward so far are now enlarged in the following 
set of experiments. 

6.5 Experiment 25. The Differentiator 

Use the three differentiating circuits of Fig. 6.8 to produce the 
derivative of the function e“ f/2 . 


c t 



The three types of differentiator shown consist of pure differentia¬ 
tion ( a ), and two approximate circuits using RC networks ( b ) and 
(c). 

1. The circuit to produce the function e 1/2 is obtained as in 



Fig. 6.9. The function required is regarded as the solution of the 
differential equation 

— + iz = 0, (6.35) 

dt 

with z = 1 at t = 0. 

Potentiometer 1 2 

Setting 0*5 1*0 
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2. In Fig. 6.8(a) the operation of the circuit is 

e 0 = -RC^p (6.36) 

dt 

and its transfer function, 

= -RCs. (6.37) 

From Eq. (6.36) it will be seen that on replacing s by jco, 

G(co) = | Y(jco)\ 

= RCco (6.38) 

and 

4>(w) = arg Y(jco) 

= -90°. (6.39) 

The gain M in decibels (db) is defined as 

M(m) = 20 log 10 G(co). (6.40) 

Substituting for G and evaluating M at co and 2co gives 


M(2m) — M(oj) = 20 log 10 2o> — 20 log 10 co 
= 20 log, o 2 
= 6 db approximately. 

This is true for all angular frequency co, and hence the gain-frequency 
plot for the differentiator increases at 6 db per octave. For example 
if RC = 1 and co = 1 then M = 20 log 1=0 db, and when co = 2 
M = 6 db. 

The phase characteristics shows a constant phase difference, the 
output lagging the input by 90 degrees. 

Unwanted noise signals are always present in analogue computers 
and they may well have a high frequency content. Because of the 
linear dependence of gain on frequency these signals will be amplified 
by a factor co when they emerge from a pure differentiator. Thus 
while their amplitude forms an insignificant part of the true signal 
on the input side this may not be so at the output. The higher the 
frequency of the noise signal the larger will be the magnitude of the 
undesired signal, thereby causing distortion of the true output. 

Another difficulty encountered with this ideal circuit is instability 
caused by using too high a capacitance in the input arm. Especially 
with computers using plug-in components, the operator should check 
that the value of C is kept to within the value stated in the manu¬ 
facturer’s rating of the d.c. amplifier. 

B.A.C.T.— 7 
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Choose RC = 1 and perform the operation y x = —(c t/2 ) and 

dt 

record y x . 

3. In Fig. 6.8(6) the transfer function is 


Y(s ) = - 


RCs 

1 + R x Cs 


(6.41) 


At low frequencies this approximates to —RCs signifying perfect 
differentiation (see §6.3 Eqs. (6.30) and (6.31)) while at high fre¬ 
quencies it becomes R(R l9 a constant gain factor. 

Choose (i) RC = 1, (ii) R X C = 0 01. 

Basically the time constant R X C must be chosen to be small 
compared with the other time constants occurring in the simulation. 
In this case R X C RC. Usually the time constant is optimized so 
as to reduce the spurious noise level to a minimum while not 
interfering significantly with the true differentiation at low fre¬ 
quencies. 

A disadvantage of this method is the constant gain at high 
frequencies. 

Perform the operation y 2 =-—— (e~ t/2 ) and record y 2 * 

1 + 0015 


4. In Fig. 6.8(c) the transfer function is 


Y(s) = 


-RCs 

(1 4* R x Cs)(l + RC x s) 


This circuit exhibits the advantages of (6) and also a diminishing 
response at high frequencies when 


The gain-frequency characteristics of the three circuits are shown in 
Fig. 6.10 for comparison, the values of the constants being taken as 

RC = 1, R X C = RC X = 001. 

The actual values of R, R u C and C x will depend on whether the 
computer used is of the fixed component or plug-in type. Care must 
be taken to ensure that any capacitor placed in series with the input 
of an amplifier is kept to a minimum, or instability may result. 
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6.6 Experiment 26. The Simple Lag Transfer Function 

Examine the frequency response of the transfer function 1/(1 + Ts) 
where T = 01 s. 

1. The circuit to produce a transfer function 1/(1 + (Ms) is 
shown in Fig. 6.11. With this arrangement, 


Y(s) = 


R 1 
/?/(! + RCs) 


(6.42) 


Choose R, R it C so that R/Ri = 1 and RC = 0-1. 

2. Calculation of Frequency Response. The steady state response 



of the transfer function to sinusoidal inputs, that is, the frequency 
response, consists of the graphs of gain and phase plotted against 
frequency (the Bode plot). 

If Cj = A sin cot then 

G(co) = IfoCM 
I ^(jco) 
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= 1 1 - 1 

1 + 0-1/co | V( 1 + O-Olco 2 ) 

and the gain in decibels is therefore 

M = —10 log 10 (1 + O Olco 2 ). 

The phase angle 

'' ,< “ ) = are (rrW) 

= —tan -1 OTco. 

For low frequencies 0*01 co 2 can be neglected compared with unity 
and then G = 1 and M = 0 db. For high frequencies OOlco 2 is large 
compared with unity and then G = 10/co and M = 20 — 20 log 10 co 
db approximately. The latter gives a straight line graph if co is 
plotted on a logarithmic scale passing through the point M = 20, 
co = 1 and has a slope of about — 6 db per octave. 

Calculate M and (j) over the range co = 1-500 rad/s and plot M 
against co on graph paper using a linear scale for M and a logarithmic 
scale for co. A useful spacing of points is obtained by evaluating 
M at co = 1, 3, 7, 10, 30, 70, 100, 300 and 500 rad/s. 

The lines M = 0 and M = 20 — 20 log 10 co are asymptotes to 
the graph. Their point of intersection is called the break point 
and the frequency at which this occurs, co r , the break frequency. 
In general this occurs when co = l/T. For the case in question, 
co c — 10 rad/s. 

Simple calculation will show that the true curve is 3 db below the 
break point at the break frequency and that, when co = \w c = 5 
rad/s, M = —1 db approx., and, when co — 2co c = 20 rad/s, 
M — — 7 db approx. These values together with the two asymptotes 
allow for quick drawing of the Bode plot for any simple lag. 

No such simple process is available for the phase diagram 
which is plotted linearly in degrees against co on the same logarithmic 
scale as in the gain plot. It will however be noticed from Eq. (6.45) 
that (j) = 0 for low frequencies, is —90° for high frequencies and is 
— 45° at the break frequency (see Fig. 6.12). 

3. Measurement of Frequency Response. Apply sinusoidal inputs 
to the lag network over the same frequency range as given above and 
record the output using a transfer function analyser. 

This unit provides the sinusoidal signal, called the reference signal. 


(6.43) 

(6.44) 

(6.45) 
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and can measure the ‘in phase’ and ‘quadrature’ components p and 
q respectively, of the output of the network under test. Tabulate p and 
q against w . 

If a is the amplitude of the reference signal (i.e. e { = a sin cot) then, 

q _ V(P 2 + q 2 ) 

a 

and 

from which the Bode plot may be constructed. 


(6.46) 

(6.47) 



a sin <ot e 0 = p sin wt + q cos wt 


Fig. 6.13 


4. Compare the results for the calculated and measured responses. 
They should agree closely. 

5. If a transfer function analyser is not available, students may use 
a low frequency oscillator to provide the sinusoidal excitation or 
indeed an oscillatory circuit using computer components (see Experi¬ 
ment 9). The adjacent channels on a pen-recorder provide a means 
of recording and comparing in amplitude and phase the input and 
output signals, a separate computer run being required for each value 
of co. 
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6.7 Experiment 27. The Quadratic Transfer Function 


Set up a circuit to produce the transfer function 20/(20 4- 4-2s + s 2 ) 
and produce calculated and measured Bode plots. 

1. The quadratic transfer function occurs frequently in the simula¬ 
tion of linear systems, e.g. the approximate response of a rate gyro. 
Its general form is 


Y(s) = 


CO 


2 

n 


co 2 + 2£c o n s + s 2 


(6.48) 


This describes how e 0 /e ; varies with frequency when e 0 and c', are 
defined by the equation 



+ 2Cco„ 


de 0 

dt 


+ co 2 n e o 


= (ofe 


(6.49) 


with £>0 = 0 = de 0 jdt at t = 0. 

2. The gain of the transfer function (6.48) in decibels is 

M(co) = 20 log jo I-- 
\e ; 


20 log t 


and the phase angle 


(j)(co) = —tan 


a /[(^« “ c ° 2 ) 2 + 

1 2£a) n Q) 


Note that for co a> n Eqs. (6.50) and (6.51) reduce to 


and 


When co > co n 


and 


M(cu) - 0 

cj>(co) = —tan” 1 « 0. 

M(co) = 40 log,o ^ 

0(co) = —tan -1 


(6.50) 

(6.51) 

(6.52) 

(6.53) 

(6.54) 


— —180° approximately. (6.55) 

Relations (6.52) and (6.54) give the low frequency and high fre¬ 
quency asymptotes to the M-co curve, the latter being a line decreas- 
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ing at 12 db per octave, while Eqs. (6.53) and (6.55) provide the 
initial and terminal values of (j>. 

For values of co near co n , the gain curve exhibits a resonance peak 
when the value of £ is small. The magnitude of this peak tends to 
infinity at co — co n when £ tends to zero. In addition, it is seen from 
Eq. (6.51) that the value of cj> at co = co n is —90° except when £ = 0. 

The foregoing discussion enables some of the characteristics of the 
Bode plot to be quickly obtained and used at least to check the 
computed curves. 

In the present example col — 20 and 2 Cco n — 4-2 giving co n = 4-47 
rad/s and £ = 0*47. Put these values into relations (6.50) and (6.51) 
and calculate M and cj) against co over the range co = 1-500 rad/s. 
Plot the curves on linear-log graph paper. 

3. By means of the method described in the previous experiment 
determine the Bode plot experimentally. The circuit of Fig. 6.14 may 
be used for this purpose. 


*3 



This circuit forms a good example of the use of complex imped¬ 
ances in nodal analysis. Using KirchhofTs first law and assuming 
an earth potential at the amplifier grid G , 


and 


ft _ fb £o_ fb 

Rl *3 


*2 


— C t se b — 0 


(6.56) 


^ + C 2 se 0 = 0 (6.57) 

R 2 


where e 0 , e h e b are the Laplace transforms of the voltages e 0 , e t and e b . 
Solving for e 0 in terms of e, gives 


gp(s) 


^3_I_ 

Rl 1 + + R 3 + -^A c 2 s + C 1 C 2 f? 2 i? 3 s 2 


(6.58) 
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If for convenience R 3 = R l then 

*ofr) = __1____ (6 59) 

<?;(s) !+(/?!+ 2R 2 )C 2 s + C 1 C 2 R 1 K 2 s 2 ' V ' 

This is to be compared with 20/(20 + 4-2 j + s 2 ) giving the relations 
OR, + 2 R 2 )C 2 = 0-21 (6.60) 

and 

C^C 2 R X R 2 = 0*05. (6.61) 

If Ci = 0-50 /iF and C 2 = 0*01 /iF then 


and 

Solving for R u R 2> 


R x + 2R 2 = 21 
R y R 2 = 10. 


= 1M, = 10M. 


(6.62) 

(6.63) 


The component values in this example have been carefully contrived. 
In general, although extremely helpful in saving computer eouipment 
and avoiding the necessity for awkward gains, this method has 
certain disadvantages: 

(a) The flexibility is decreased as far as system parameter changes 
are concerned. 

(b ) The complex relations between the component values and 
parameters often call for computational complexity. 

0 c ) Unusual component values frequently result from ( b ). The use 
of nomographs relating circuit components to parameter 
values alleviates the second disadvantage. 

Table 6.1 gives values of M(oS) and 0(ca) from which the Bode 
plot for the given transfer function may be obtained. 


6.8 Experiment 28. The Synthesis of Transfer Functions 

In a certain servo-mechanism the transfer function between the 
Laplace transforms of the output variable 6 0 and the error signal e is 


9 0 = 0*5 

£ s(l -T 0*2^)(1 + 5 ) 


(6.64) 


Set up this transfer function using the RC network method and 
test it for frequency response. 

1. It has already been mentioned that linear transfer functions 
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TABLE 6.1 


co (rad/s) 

1 

3 

7 

10 30 

M(co) 

0-238 

1-553 

-6-297 

13-098 -32-958 

<t> M 

-0-2176 

-0-8532 

-2-3495 

-2-6581 -2-9993 

co (rad/s) 

70 

100 

300 

500 

M(a>) 

— 47*764 

-53*969 

-73*062 

-81-937 

Hu ) 

-3-0811 

-3-0995 

-3*1276 

-3-1332 


may be multiplied together and manipulated algebraically. Thus 
relation (6.64) is equivalent to four separate blocks as shown in 
Fig. 6.15. 


i i i 



12 3 4 


Fig. 6.15 


The response to sinusoidal inputs is as usual obtained by replacing 
s by jco in Eq. (6.64) and similarly in each separate block of Fig. 6.15. 
Table 6.2 shows G(co) and for each block. 

The term 0*5 has gain in decibels 20 log 10 0-5 = —6 and is con¬ 
stant. The phase angle is also constant at zero. The graph corres- 


TABLE 6.2 


Block 

1 

2 

3 

4 

G(a>) 

0*5 

1 

1 

1 


co 

v'(l + 0-04CO 2 ) 

V(1 + g* 2 ) 




— tan 1 co 


0 


-90 


— tan 1 0*2 co 
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ponding to the term 1 /s has a gain slope of — 6 db per octave and 
cuts the zero db line when co = 1. The phase angle is constant at 
“90°. The graphs corresponding to both the terms 1/(1 + s) and 
1/(1 + 0*2^) have M = 0 db as low frequency asymptote and have 
a high frequency asymptote with slope — 6 db per octave. The former 
has a break point at co = 1 rad/s and the latter at co = 5 rad/s. At 
these frequencies the true response curve lies 3 db below the inter¬ 
section of the asymptotes. 

Figure 6.16 shows the Bode gain plot for each block together with 
the overall response obtained by adding together the individual M 
components. The curves are labelled according to the numbering of 
the blocks in Table 6.2. The overall Bode phase plot may be obtained 
by adding together </> values for each block in a similar way. 


Frequency w(rad/s) 



Fig. 6.16 


2. To simulate the transfer function the following circuit can be 
used. 



The values of resistance and capacitance are shown in megohms 
and microfarads respectively. In fixed gain computers it may be 
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necessary to utilize different values of these components, when the 
following relations must be satisfied. 

R 1 C 1 = 2, \ 

R 3 /R 2 = 1 = R 5 IR 47 (6.65) 

R ?>^2 = 0 * 2 , | 

R 5 C 3 = 1. j 

For instance in a fixed component transistorized computer the value 
of the capacitor may be 10 fj. F. 

Then C\ = C 2 — C 3 = 10, and it follows that R t = 0*2, 
r 2 = 0 02 = R 3 and R 4 = 01 = R 5 . 

3. Perform the frequency response test to determine the Bode plot 
for 6 0 ls and check the graphs by using theoretical values. 

6.9 Experiment 29. Simulation of Automatic Control 

Set up a circuit to investigate the effects of step changes in the desired 
output variable of a control system having a process represented by 
the transfer function 

= Y(s ) =- - - ( 6 . 66 ) 

K(s) (1 + 01s)(l + 0-25s) 

and governed by a three-term controller. 



Fig. 6.18 


1. The object of an automatic control system is to hold a measured 
variable such as temperature or pressure at some desired value. 

This is achieved by comparing the actual value, 0 o , with the 
desired value, and producing appropriate action, F, from the 
obtained error signal 

e e = 0 o - 0;. (6.67) 

For example, if in a chemical process the desired temperature is 
100° C and the measured temperature falls to 95° C, a controller 
acts in such a way that more heat is provided and the temperature 
increases. If the measured temperature rises above 100° C the 
controller reverses its action. 
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2. A controller can have several physical forms, e.g. pneumatic, 
hydraulic, electrical, depending in kind on the nature of the con¬ 
trolled variables. 

There are three types of control action: (a) Proportional, (b) In¬ 
tegral, (c) Derivative. 

In (a) the action of the controller V is proportional to the error 
signal i.e. 

V = K P 0 S . (6.68) 

In (b) the rate of action, dV/dt, is proportional to 0 E , i.e. 


or 


dV 

dt 


KA 


(6.69) 


V = Kj 


e F dt . 


(6.70) 


In (c) action is proportional to rate of change of error signal, i.e. 


V = 



(6.71) 


Transforming Eqs. (6.68), (6.69), and (6.71) and assuming zero 
initial conditions, the complete action of the controller in transfer 
function form is 


V = 



+ K D .S'j(J £ . 


(6.72) 


This is shown schematically in Fig. 6.19. 





Fig. 6.19 


3. Using Eqs. (6.67), (6.68) and (6.72) the following circuit may be 
obtained to represent the control system. 
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Amplifier number 1 is used as a differencing element to simulate 
Eq. (6.67), amplifiers 2, 3 and 4 constitute the controller with an 
extra gain of 10 provided in the summing amplifier number 5, and 
finally amplifiers 6 and 7 represent the process dynamics. 



It will be seen that an approximate form of differentiation is used 
(amplifier 4). The transfer function of this has already been discussed, 
as have the simple lag circuits of amplifiers 6 and 7. 

4. Representative values for the three controller constants are 

0 ^ K P S 20, \ 

OS Kj S 100, (6.73) 

0 ^ K d S 25. J 

Thus potentiometers 1, 2 and 3 are set to record A^/lOO, KpjlO and 
K D j25 respectively. A gain of 10 is incorporated in amplifier 5 
leaving amplifiers 2, 3 and 4 with gains of 10, 2 and 2-5 respectively. 

In amplifier 4 this amounts to choosing values of R i and C t such 
that 

R.C, = 2*5 
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and ensuring that RC t is small in comparison with the lowest time 
constant, 0-1 s. A value of about 0 001 may be tried. 

Further, choose values of R 2 , R 3 , R 4 , R 5 and C 2 , C 3 to satisfy 
the following relations 

j _ ^5 

R 2 R4 

RtC 2 = 01 , 

r 5 C 3 = 0-25. 

5. Initially it is assumed that with 9 t = Oand0 o = 0 the controller 
is acting to keep the actual output equal to the desired output. Any 
input signal 0* will be regarded as representing a change in the 
desired output. Investigate and record the following modes of control 
using step changes in 0 £ to represent sudden changes in desired out¬ 
put. 

(a) No integral or derivative action ( Kj , K D both zero) using small, 
medium and large amounts of proportional action (three 
runs). Notice how as K P is increased, the steady state error 
(known as ‘offset’) is decreased whilst the output tends to 
become more oscillatory. 

(b) With each of the above use small and large amounts of integral 
action (six runs). Notice that the offset is now eliminated at the 
expense of stability. 

(c) Use small and large amounts of derivative action with each 
of the above (eighteen runs). Notice how this improves 
stability. 


\ 

(6.74) 


6.10 Further Exercises 

1. Differentiation can be performed using the implicit function 
techniques rather than by RC circuits. Simple analysis of Fig. 6.21 
will show that 

e 0 (s) __ 5 

efe) 1 + sT 

The advantages are 

(a) Only standard computing units are used. 

(b) The appropriate sign of the derivative is generated. 

(c) The break frequency l/T is easily controllable and prior 
knowledge of the noise frequencies is not necessary. 
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Use this circuit to differentiate y = sin 4 1 for various values of T. 

2. Set up a circuit to simulate the transfer function 

y( s ) =---. 

s(l + 01s) 

Obtain the Bode diagrams for this circuit experimentally and check 
them by calculation. 

3. If a function of time t,/X0> is delayed by T seconds the resulting 
function f 0 (t) is given by 

m = m - n 

It may be shown that the transfer function representing the delay is 




where s is the complex operator previously mentioned. 
Replacing s by yco for sinusoidal inputs the gain 


GfoO- 


7oO)| 


== 1 (constant) 


and the phase 


4>(co) = arg 


/oO) 


= -coT 

which is proportional to frequency. Now 


= 1 - sT + 


s 2 T 2 s 3 T 3 


which may be approximated to by the expression 
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j - 2I T 

s~+J/T 


l — sT + 1JL 


s 2 T 2 5 3 T 3 


+ 


Significant errors arise in the fourth and succeeding terms if sTis not 
small. The following circuit (Fig. 6.22) has -(s - 2/T)/(s + 2/T) 
as transfer function. Use it to test the delay of functions 


(l)/* = sin 2/, (2) fi = t 

when T = 0-25 s. 



(S - 2/T) 
is + 2/T)' 


./* 


Circuits of more complex form can be produced to provide much 
better approximations to time delays. 

4. In the coplanar study of the homing of a guided missile steering 
under the proportional navigation system the equations may be 


Reference axis 
For x T - 


Collision Point 

♦ r 


Target 


Line of sight 



_ — -♦ Initial line of sight 
Reference axis 


Fig. 6.23 


reduced to the following simplified form 


6 = Y(s)& 
Xt~ X m 


a —- 


R 


R = R 0 - tV cos 
X M = V& cos 0 O 
X T as desired 
(0 O , R 0 are the initial values of 0 and R) 


Missile rate of turn 

Line of sight angle 

Range (missile to target) 
Missile lateral acceleration 
Target lateral acceleration 



THE USE OF TRANSFER FUNCTIONS 


101 


7(5) is a transfer function representing the missile dynamics. It 
dictates the ability of the missile to perform the programmed steering 
equation 6 = Kg. 

In this problem take 


Y(s) = 


K 

1 + Ts 


where 0 < K < 10 and T = 0-5. 

Examine the miss distance (the value of X T - X M when R = 0) 
for values of K equal to 01, 1, 2 and 10 for a 2g target turn, i.e. 
X f = 64 ft/s 2 . 

The maximum values of the variables are given in the Table 6.3 
and a suggested circuit in Fig. 6.24. 

The experiment should be performed for various initial heading 
errors, i.e. ( X M - X T ) 0 . 

If a constant signal is applied through an extra input to amplifier 1, 


TABLE 6.3 


e 

& 

G 

R 

R 

1 rad 

10 rad/s 

0T rad 

20,000 ft 

1000 ft/s 

Xm 

Xm 

X T 

X T 

Xt ~ X M 

2000 ft 

1000 ft/s 

2000 ft 

1000 ft/s 

1000 ft 


Potentiometer 

1 

2 3 

4 


Setting 

KI10 1/10T Kcos0 ° 

M X T ) Q 



1000 

1000 


Potentiometer 

5 

6 7 



Setting 

X T 

V cos 0 o R 0 




1000 

20,000 20,000 



b.a.c.t.—8 
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the effect of a bias in the measurement of sight line angle may be 
investigated. 

Amplifiers 1 to 3 represent the transfer function 0*ls/(l + Ts ) 
in feedback form. A more realistic representation of the missile 
dynamics is obtained by the inclusion in the circuit of a term 
1/[1 + (2£s/c0„) + 0/co„) 2 ] immediately prior to potentiometer 1. 
The use of RC networks to simulate these transfer functions is of 
course acceptable. 
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